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A bosonic analogue of the fractional quantum Hall effect occurs in rapidly rotating trapped 
Bose gases: There is a transition from uncorrelated Hartree states to strongly correlated states 
such as the Laughlin wave function. This physics may be described by effective Hamiltonians 
with delta interactions acting on a bosonic iV-body Bargmann space of analytic functions. In 
a previous paper [RSY] we studied the case of a quadratic plus quartic trapping potential and 
derived conditions on the parameters of the model for its ground state to be asymptotically 
strongly correlated. This relied essentially on energy upper bounds using quantum Hall trial 
states, incorporating the correlations of the Bose-Laughlin state in addition to a multiply 
quantized vortex pinned at the origin. In this paper we investigate in more details the density 
of these trial states, thereby substantiating further the physical picture described in [RSY], 
improving our energy estimates and allowing to consider more general trapping potentials. Our 
analysis is based on the interpretation of the densities of quantum Hall trial states as Gibbs 
measures of classical 2D Coulomb gases (plasma analogy). New estimates on the mean-field 
limit of such systems are presented. 
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1 Introduction 

The advent in the 90 's of powerful techniques to cool and trap atoms has opened the way to 
new investigations of quantum phenomena on a macroscopic scale. It has now become possible to 
isolate extremely cold and dilute atomic gases and maintain them in mctastable states that can 
be efficiently modeled as ground states of effective many-body Hamiltonians with repulsive short 
range interactions. Although the atoms typically trapped are neutral, one can impose an artificial 
magnetic field to them. This is usually achieved by rotating the trap |Co[ IFetj although there now 
exist more refined techniques |DGJOj . Exploiting the analogy between the Lorentz and Coriolis 
force, one can easily realize that the Hamiltonian for neutral atoms in a rotating frame resembles 
that of charged particles in a uniform magnetic field. 

The analogy has been demonstrated by the nucleation of quantized vortices in cold rotating 
Bose gases (see [Fetj and references therein, in particular |BSSDj ). similar to those appearing 
in type II superconductors submitted to external magnetic fields. In this regime, the atoms all 
condense in the same one-particle state and the gas forms a Bose-Einstein condensate (BEC). 
The peculiar properties of quantized vortices, demonstrating the supcrfiuid nature of BECs, have 
motivated numerous theoretical and mathematical works, see |Aftl IFetl ICol ICPRY11 ICPRY2] and 
references therein. 

An even more striking possibility is to create with cold atomic gases phases characteristic 
of the fractional quantum Hall effect (FQHE) , a phenomenon originally observed in 2D electron 
gases submitted to very large magnetic fields jSTGj . In this regime, the Bose gas is no longer a 
condensate and mean-field theories fail to capture the physics of the system. One has to use a truly 
many-body description, and strongly correlated phases, such as the celebrated Laughlin state |Lauj 
may occur. This regime has been hitherto elusive in trapped Bose gases. For large rotation speeds, 
the centrifugal force can compensate the trapping force and lead to an instability of the gas. As of 
now, the competing demands of a large rotation speed on the one hand, necessary for entering the 
FQHE regime, and of a sufficiently stable system on the other hand, have left the FQHE regime 
unattainable with current technology (see |RRD) for a precise discussion of this point). 

Several ideas to get around this difficulty have been proposed, and it is one of those that we 
examine from a mathematical point of view, here and in the companion paper [RSY] . As we will see, 
the proposed modification of the experimental set-up also leads to new physics and mathematics. 

When discussing the quantum Hall regime it is common jAftl IFetl ICoj to assume that the gas 
is essentially 2D, that the single-particle states are restricted to lie in the ground eigenspace of the 
"magnetic" kinetic energy operator (i.e. in the lowest Landau level) and that the inter-particles 
interactions are described by a Dirac delta potential. We will give more details on these approx- 
imations, that can to some extent be backed with rigorous mathematics |LeSej . in Section 12.11 
Assuming their validity, the system may be described in the rotating frame by the following many- 



Sj N = \ F{z u . . . ,ZAr)e" E f=i |Zj|2/2 € L 2 (R 2N ), F holomorphic and symmetric } . (1.2) 



By symmetric we mean invariant under the exchange of two particles Zi and Zj for i ^ j, i.e. 
bosonic. Here O > is the rotation speed, with the convention that the rotation vector points 



body Hamiltonian : 




(1.1) 



acting on the Lowest Landau Level (LLL) for N bosons 



1 h stands for Landau, referring to the lowest Landau level. 
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in the direction perpendicular to the plane in which the particles are confined. The interaction 
between particles i and j is given by the delta function gS(xi — Xj) where g is proportional to the 
scattering length of the original 3D interaction potential between the particles fT.iSe . In our units 
H = m = 1 and we will also fix O = 1 in the sequel. 

To find the ground energy state E h of the system we minimize the expectation value (iJ L )^ 
amongst iV-particles states \1/ in Sj N : 

E L := inf |£ L [*] = <*, H L *) , * G fi N , ||¥|U»(r»w) = l} . (1.3) 

More precisely, the energy of a state ^> reads 

N r ( i 



, j . . / I V / 1 / , 

J=l K V 7 l<i<j<AT 



^[*]:=V/ ^)-o n N |*W^+<? V / mzi = z i )\'dZ i (1.4) 



2(iV-l) 



where Z = (zi, . . . , z/v), c?Z = dz\ . . . dzrf, ^(z; = Zj) is the function evaluated on the diagonal 
Zj = Zj and dZi is e?Z whith the factor dzi removed. Note that it makes perfect sense to use a 
delta interaction potential in S) N : since functions in this space arc all smooth, the integral 

S(zi - z 2 )*> := / \^{z 2 ,z 2 ,z 3 ,...,z N )\ 2 dz 2 ...dz N = / |*(zi = z 2 )| 2 dZ! 

is always well-defined. 

In most experiments, the trapping potential V is quadratic 



and hence the effective potential, 



V(x) = ±tf ± \x\ 2 



V cS {x) = V{x)-\n 2 \x\ 2 , (1.5) 



that includes the effect of the centrifugal force is bounded below only if f2 < f2j_. There is thus a 
maximum velocity that one can impose on the gas, and the FQHE regime of strong correlations is 
expected to occur when f2 — > i.e. when the centrifugal force almost compensates the trapping 
force. To avoid this singularity of the limit O — > Qj_ , it has been proposed |MFi IViej to add a weak 
anharmonic component to the trapping potential. A simple example is provided by a quadratic 
plus quartic trap of the form 

V(r) = ^n 2 r 2 + kr 4 . (1.6) 

With such a potential there is no theoretical limit to the rotation speed that one can impose on 
the system. One may thus expect the regime Q — > Cl± to be less singular and more manageable 
experimentally. This has been demonstrated in the Gross-Pitacvskii regime by using this kind of 
trap to rotate Bose-Einstein condensates beyond the centrifugal limit [BSSDj . 

In |RSY] and the present paper we focus on the particular case where the trap V in ([l.ip is 
given by (|1.6[) : 

N 

f l = ]T H^l 2 + fcN 4 ) +g y Zs(z i - zj) (1.7) 

3=1 i<3 

with 

ui-.= l (n 2 ± - n 2 ) . 
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We are mainly interested in two questions : 



4 



1. In which parameter regime can one obtain strongly correlated states as approximate ground 
states for this model ? 

2. Can the addition of the quartic part of the potential lead to new physics and to phases not 
accessible with the purely harmonic trap ? 

Let us now be more precise about what we mean by strongly correlated states. In the simpler case 
where k = 0, it is easy to see (we give more details in Section HO]) that for sufficiently large ratio 
g/uj, the exact ground state of the system is given by the bosonic Laughlin state 

* Lau = CLau II {Zi-ztfe-^^l* (i. 8) 

l<i<j<N 

where CLau is a normalization constant. This wave function^] was originally introduced in |Laul 
ILau2j as a proposal to approximate the ground state of a 2D electron gas in a strong magnetic 
field. The correlations encoded in the holomorphic factor in (|1.8[) decrease the interaction energy 
and even cancel it in our case of contact interactions. 

More generally (see |LeSel Section 2.2] for further discussion), it is expected that for g/{Nuj) 
of order unity, one encounters a series of strongly correlated states with smaller and smaller in- 
teraction energy, ultimately leading to the Laughlin state. Several candidates have been proposed 
in the literature (see [Col IViej for review) for the other correlated states that should occur before 
eventually reaching the Laughlin state. Although there is some numerical evidence that they have 
a good overlap with the true ground states, none of these trial functions is as firmly established as 
the Laughlin state. 

On the other hand, for small g/ (Nui) one can show that the ground state energy of (|1.7p is well 
approximated by Gross-Pitaevskii theory |LSYj . that is by taking a Hartree trial state 

N 

*( Z1 ,..., ZN )= n^fe) (°) 

for some single particle wave function <f>. The corresponding GP theory with states restricted to 
the lowest Landau level has some very interesting features, studied in jABNli [ABN2j . One thus 
goes from a fully non correlated state for small g/(Nui) to a highly correlated Laughlin state for 
large g/{Nu). 

Now, what changes in this scenario when a quartic component is added to the trap? The tools 
of |LSY| still apply in the Gross-Pitaevskii regime |Draj . which leads to a GP theory with some 
specific new aspects |BR[ IR] . The strongly correlated regime is more difficult: The Laughlin state 
is no longer an exact eigenstate of the Hamiltonian, so it is not obvious that one should eventually 
reach this state. One should certainly expect that the ground state will essentially live in the kernel 
of the interaction operator 2jy = J2i<j $( z i ~ z j) 

Kev(l N ) = {f Lau F(zi, . . . , z N ), F holomorphic and symmetric} C fi N (1. 10) 

in a regime where the interaction energy dominates the physics, i.e. when g is sufficiently large 
and uj, k sufficiently small. It is of interest to derive explicit conditions on the order of magnitude 
of these parameters that is needed to obtain a strongly correlated state in a certain limit, and also 
to obtain rigorous estimates for the energy E h in this regime. 

Also, a remarkable new feature of the model with k =/= is that it is now possible to consider 
negative values of u), that is to go beyond the centrifugal limit of the harmonic trap by taking 

2 or rather its fermionic analogue where the exponent of the (z^ — zj) factor is an odd number, ensuring fermionic 
symmetry 
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> Qj_. In this case, the effective potential (|1.5[) has a local maximum at the origin, which can 
lead to new physics. Indeed, it can be shown (see below) that the matter density of the Laughlin 
state is almost flat in a disc around the origin. If the trap develops a maximum there, the Laughlin 
state will have a large potential energy and one may expect that other strongly correlated phases 
belonging to Ker(Ijv), with depleted density at the center of the trap, will be preferred to the 
Laughlin state, even for large g and small uj, k. 

Our main results on the model (|1.7p . derived in |RS Yj . can be summarized as follows 

• we rigorously identify conditions on the parameters of the problem g, uj, k and N under which 
the ground state of H h is fully correlated in the sense that its projection on the orthogonal 
of (|1.10p vanishes in a certain limit. 

• within this regime we obtain upper and lower bounds to the energy whose order of magnitudes 
match in the appropriate limit. 

• we can identify a regime where the Laughlin state does not approximate the true ground 
state, even if it almost fully lives in Ker(ZAr). We prove that a state with higher angular 
momentum is preferred, which can never happen in a purely harmonic trap with k = 0. 

Precise statements are given in Section [2] below, after we have recalled several facts about the 
formulation of Problem (|1.3I) using the Bargmann space of holomorphic functions. The proofs are 
given in |h'SV| Section 5]. They mostly rely on adequate energy upper bounds derived with fully- 
correlated states of the form (|1.10[) . Indeed, since the Laughlin state is no longer a true cigenstate 
of the Hamiltonian, the evaluation of its potential energy, and that of other candidate trial states, 
is a non trivial problem. Also, as discussed above, the effective potential changes from having a 
local minimum at the origin to having a local maximum when uj is decreased. We thus need to 
have some flexibility in the matter density of our trial states to adapt to this behavior, which leads 
us to the form 

*£ = c m f[z? J] {z,-z 3 fer^U\^/- (1.11) 

j = l l<i<j<N 

were c m is a normalization constant. One recovers the pure Laughlin state for m = and the 
states with m > are usually referred to as Laughlin quasi- holes jLaui ILau2] (hence the label qh). 
They were introduced as elementary excitations of the pure Laughlin state. The factor rX,=i Z T * s 
interpreted as an additional multiply quantized vortex located at the origin. Its role is to deplete 
the density of the state when uj < to reduce potential energy. 

We have no rigorous argument allowing to prove that the true ground states are asymptotically 
of the form above in the strongly correlated regime although we do believe that it is the case. 
A rudimentary lower bound to the energy confirms however that the trial states (|1.11[) at least 
give the correct order of magnitude for the energy in the strongly correlated regime, as stated in 
Theorem 12.31 below. 

The simple method we use in |RSY] to evaluate the energy of these states only allows a limited 
physical interpretation of the results. It is also not very flexible and docs not apply to a very large 
class of trapping potentials. Indeed, it relies on the representation of the effective potential 

V u , k (r) = ujr 2 + kr A (1.12) 

using angular momentum operators on the Bargmann space of analytic functions (see Section [2.11 
below). This requires that the potential be radial and hence does not allow to treat anisotropic 
potentials. Also, the estimates detailed in jRSYl Section 5] rely on the fact that the quartic term 
in the potential comes with a positive coefficient in front. Even for a simple potential like 

V(r) = -k 4 r 4 + k e r 6 
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with k&, ke > 0, the method fails. We shall not treat the minimization problem (|1 .3[) with such a 
potential because using the method of the present paper it would be straightforward generalization 
of the case we treat, only with more tedious computations. In any case, we consider it important 
to have a method that works for more complicated potentials than (|1.12[) . 

Since we are interested in the strongly correlated regime, all the trial states we use will belong 
to Ker(ZAr) and thus have zero interaction energy. Our task is then to calculate their potential 
energy 

N [ V{z)p^(z)dz (1.13) 
Jm 2 

where 



py(z)=l \^{z,z 2 ,...,z N )\ 2 dz 2 ...dz N (1-14) 

J R 2(JV-1) 

is the one-particle density of the state VP, that we have defined so that its integral is 1. Given a 
candidate trial state we thus want to evaluate precisely what the corresponding matter density is. 
This will also provide a better understanding of the wave functions p. Ill) that play a central role in 
FQHE physics [Girl ILFSl ISTGj . As we will see, this understanding is crucial for the interpretation 
of our results on Problem 11.31 

Our main tool is the well-known plasma analogy, originating in ILau|lLau2] . wherein the density 
of the Laughlin state is interpreted as the Gibbs measure of a classical 2D Coulomb gas (one 
component plasma). More precisely, after a scaling of space variables, one can identify the N- 
particle density of the Laughlin state with the Gibbs measure of a 2D jellium with mean-field 
scaling, that is a system of N particles in the plane interacting via weak (with a prefactor iV -1 ) 
logarithmic potentials and with a constant neutralizing background. Within this analogy the vortex 
of degree m in p. lip is interpreted as an additional point charge pinned at the origin. 

Existing knowledge (e.g. jCLMPl iKiel iKSl iMSl 155] ) about the mean- field limit for classical 
particles then suggests that we should be able to extract valuable information on the density of 
the Laughlin and related states in the limit N — > oo. For our purpose, precise estimates that are 
not available in the literature are required, so we develop a new strategy for the study of the mean 
field limit that gives explicit and quantitative estimates on the fluctuations from the mean field 
density that arc of independent interest. 

It is convenient to work with scaled variables, defining (we do not emphasize the dependence 
on m) 

p N (Z) :=N N \*£(VNZ)\ 2 . (1.15) 
The plasma analogy consists in comparing the one-body density corresponding to (|1.15[) , i.e. 

(J$(z)= / fi N (z,z 2 ,...,z N )dz 2 ...dz Nl 

JR2(«-1) 

with the minimizer g MF of the mean-field free energy functional 

£ MF [p]= [ W mP + 2D(p,p) + N- 1 [ plogp (1.16) 

JR 2 JR 2 

amongst probability measures on R 2 , p £ P(R 2 ). Here 

W m (r) = |rf-2^1ogr 

and the notation 

D(p,p) = - p{x)\og\x-y\p(y)dxdy (1.17) 
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stands for the 2D Coulomb energy. 

We will provide new estimates for this classical problem, discussed at length in Section [3] 
Theorem 13.21 therein is our main result in this direction. Combined with Proposition 13. II below it 
can be summarized as follows 

Theorem 1.1 (Plasma analogy for QH trial states). 

There exists a constant C > such that for large enough N and any smooth function V on M 2 



V (j$ - g MF ) < CN- 1 ' 2 log N\\VV\\ L 2 {R2) + CN^W] 



if 'to < N 2 , and 



ifm^>N 2 



KCN-^m-^WVW 



(1.18) 



(1.19) 



This result allows to evaluate the one-particle density of our trial states and deduce estimates 
of the potential energy (|1.13[) . Optimizing these over to, we find 



'-opt 



if u > -2kN 

-jfc -N if u < -2kN. 



(1.20) 



We interpret this as a strong indication that, within the fully correlated regime, a transition 
occurs for lo < and |w| oc kN between a pure Laughlin state and a correlated state with a density 
depletion at the origin. Interestingly we also find that the character of the mean-field density g MF 
and thus that of the one-particle density of the state (jl.lip strongly depends on to: For m -C N 2 it 
is correctly approximated by a flat density profile located in a disc or an annulus (depending on the 
value of to), whereas for to 3> ./V 2 , the density profile is approximately a radial Gaussian centered 
on some circle. This is due to a transition from a dominantly electrostatic to a dominantly thermal 
behavior of the 2D Coulomb gas to which we compare our trial states. Using the expression of the 
optimal value of m given in (|1.20[) this suggests a further transition in the ground state of (|1.3p 
in the regime \oj\ oc kN 2 . Establishing these phenomena rigorously remains a challenging open 
problem. 



The rest of the paper is organized as follows: In Section 12.11 we formulate (|1.3[) precisely as 
a minimization problem over the Bargmann space of analytic functions. In Section 12.21 we state 
our main results about this model, whose proofs are given in |RSYj , and discuss in Section 12.31 
how Theorem 11.11 allows to interpret and improve them. Section [3] is the mathematical core of the 
present paper. It contains the details on the plasma analogy, along with the proof of Theorem ll.il 
and related results. Finally, Section [4] shows how to use the plasma analogy to obtain rigorous 
estimates of the potential energy of our QH trial states. 



2 Rotating bosons in the lowest Landau level 
2.1 The model 

In this section we discuss the derivation of the effective Hamiltonian that we will study and its 
most important properties. Our starting point, Equation (|2.1[) below, is the full 3D Hamiltonian 
for a rotating trapped Bose gas with repulsive interactions given by a two-body potential v > 0. 
We make three standard approximations: (1) the motion along the axis of rotation is frozen, (2) 
the states for the motion in the plane are reduced to the lowest Landau level, (3) the interaction 
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potential is replaced by a delta function gS where g is proportional to the scattering length of v 
(see [LSSY1 Appendix C] for a definition). 

These approximations have been studied in [LeSej in the case k = 0. Rigorous bounds quantify- 
ing their validity have been derived, and the approach can be adapted with almost no modifications 
to the case k ^ 0. For these reasons we only sketch the derivation of the reduced Hamiltonian and 
focus on its essential properties. 

The full many-body Hamiltonian 

The rotating Bose gases where it is proposed to try to create quantum Hall phases can be described 
in the rotating frame using the following many-body Hamiltonian 

N n ^ ^ 

i i ^ ' • j 

Here Xj = {xj,i,Zj,2,Zj,3) G K 3 is the coordinate of the j-th particle, V is a confining external 
potential and v the two-body interaction potential. Units are chosen so that ft = m = 1. We also 
choose the coordinate axis so that the rotation vector £1 = ^le^ is proportional to the unit vector 
e3 in the 3-dircction. The vector potential 

A(x) =Q(x 2 ,-x 1 ,0) 

in (|2.1[) represents the Coriolis force (analogous to the Lorentz force) while the negative quadratic 
potential proportional to fl 2 corresponds to to the centrifugal force. We are interested in the 
ground state of this Hamiltonian. 

We consider a trapping potential of the form 

V(x) = ifii (x\ +xj)+k (xj + xlf + |n N a^ (2.2) 

and it is necessary for Hff to be bounded below to require that the effective potential 

K ff (x) = V(x)-±n 2 (x 2 1+ x 2 2 ), (2.3) 

taking into account the effect of the centrifugal force, remains bounded below. This is ensured if 
either k > or k = and f2 < 

Reduction to the lowest Landau level 

The 1-particle Hamiltonian in (|2.1[) is given by 

H 1 = \ (iV + A(x)) 2 + \{-d 2 + nf x 2 ) + y(x) - ^(x 2 + x\). (2.4) 

The first term is the 2D Landau Hamiltonian with spectrum 2(n + i)0, n e N, the second a 
one-dimensional harmonic oscillator in the 3-direction with spectrum (ny + 5)^11 , Till € N. These 
two terms commute and can thus be diagonalized simultaneously. 

When the energy scales associated with the effective trapping potential in the 12-plane 
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are much smaller than the gaps between the energy levels of the first two operators in (j2.4[) it is 
natural to restrict attention to joint eigenstates of the first two terms of (|2.4p with n = n\\ = 0. 
The motion in the 3-direction is then 'frozen' in the ground state of the harmonic oscillatoo As 
far as the 12-plane is concerned the state is in the lowest Landau level (LLL). 

Henceforth we choose units so that fi = 1. Replacing (xi,X2) by the complex variables z = 
Xi + 1x2 and z = x\ — 1x2 and denoting 8 Z = \{8\ — 182), 8% = 8 Z = \{d\ + 182) we can write the 
Landau Hamiltonian as 

i(iV ± +A(x)) 2 =2(V a + i) (2.5) 

with the creation and annihilation operators a) = —8 + ^z, a = 8 + ^z. Eigcnfunctions ip(z, z) in 
the lowest Landau level are solutions of the equation aip = 0, i.e., 

8 z ip(z,z) = --zip(z,z). 

They are therefore of the form 

ifi(z,z) = /(z)exp(-|z| 2 /2) (2.6) 
with 8 z f(z) = 0, i.e., / is an analytic function of z. 



The Bargmann space 

As seen above, single particle wave functions in the LLL correspond uniquely to functions in the 
Bargmann [Bar! IGJj space B of analytic functions / on C such that 

(f,f)B= [ |/(z)| 2 exp(-|z| 2 )dz<oo (2.7) 
Jc 

where dz denotes the Lebesgue measure on C ~ M 2 . We denote by Sj the space of the full wave 
functions ip including the Gaussian factor in (|2.6|) . It is a subspacc of the Hilbert space of square 
intcgrable functions w.r.t. dz. Thus, state vectors in the LLL can either be regarded as elements 
of B or of Sj and we shall make use of both points of view. 

For our TV-body system of bosons in the LLL the corresponding Hilbert spaces are the symmetric 
tensor powers of B or fj, denoted by B N and Sj N : 

B N = {f holomorphic and symmetric such that F(zt, z N )e~ ,2 E L 2 (R 2JV )| (2.8) 

% N = ...,z N ) = F( Zl , . . . , z N )er Ef=i l^l 2 / 2 G L 2 (R 2W ), F e B N ) . (2.9) 

The scalar product on B N is given by 

(F, G) BN = (Fe-^M'/', Ge -Sr=il«il 8 /2\ . (2 . 10) 



The reduced TV-body Hamiltonian 

We can now define our energy functional. When the above reductions have been made, the only 
term in the one-particle Hamiltonian (|2.4[) that is not fixed is the effective potential term in the 12 
plane. For a short range potential it makes sense physically to replace v by a delta potential g8. It 
also makes sense mathematically since the wave functions of the lowest Landau level are smooth. 

3 We could also consider a more general trapping potential in the 3-dircction. 
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Justifying the substitution rigorously is a difficult task |LeSej on which we will not elaborate, but 
if we take this for granted we obtain 

£ L [V]=N [ V u , k (r)p*{z)dz + AN(N-l)g [ |tf(z 2) z 2 , z 3 . . . , z N )\ 2 dz 2 . . . dz N (2.11) 

JR2 J R 2(]V-1) 

where pqr is the one-body density normalized to 1, g the coupling constant and 

V u , k (r) =ur 2 + kr\ (2.12) 

The ground-state energy E L is then defined as in (jl.3p . It is useful to reformulate this problem in 
the Bargmann space: 

Lemma 2.1 (Hamiltonian in the Bargmann space). 

Define the Hamiltonian acting on B N 

N 

H B :=N(uj + 2k)+^2((uj + 3k)L ] +kL 2 J )+g ^ 5 l3 (2.13) 

j = l l<i<j<N 

with Lj = Zjd Zj the angular momentum operator in the j-th variable and 

S l3 F(. ..,z i ,...,z j ...) = —F( ...,^(zi + zj), -(Zi + Zj ), ...). (2.14) 



We h 

£ L [^] = (F,H B F) BN w ith* = Fexp^-J2^] , F€ B N . (2.15) 
In particular 

E L = E B := inf (t b nH b . (2.16) 

Proof. The fact that the 5 interaction operator acts as (|2.14[) seems to have been noticed first in 
jPB| . Indeed, taking (|2.14[) as the definition of the interaction operator and using the analyticity 
of F € B N , one easily realizes that 

(F,%F) BN = f \F(--- ,z,--- ,z,---)\ 2 eM-2\z\ 2 )dz 

x eM-J2 k ^k\ 2 )I\ k ^dz k . (2.17) 
It is sufficient to perform the computation for a two-body / S B 2 : 

</, S 12 f) B2 = ± JJ c2 f( Zl ,z 2 )f (*±*, *±*) e-M'-M'dKd* 

1 /(ii,w)e- 2|M|2 / f(u + v,u- v)e- 2 ^ 2 dvdu. 
Jc 

Fixing u, writing g(v) = f(u + v, u — v) with an holomorphic g(v) = J2n>o 9nV n we find 

f{u + v,u-v)e- 2 ^dv=g f e^dv = Jp(0) = u), 
c Jc z 1 

which is what we need. 
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On B the operator of the 3-component of the angular momentum is L z = zd z and an integration 
by parts shows (see e.g. jABl Lemma 3.1]) that, for / £ B, 



(f,L z f) B = / (\z\ A - l)\f{z)\ z cM-\zV)dz. (2.18) 
Jc 

Integrating by parts twice we also obtain 

\z\ 4 \f(z)\ 2 exp(-\z\ 2 )dz = (f,L 2 J) B + 3(f,L z f) B + 2(f,f) B . (2.19) 



Putting ([2~T7)) . ([2~T8)l and (|2~T9)l together proves ([2~T5]) and (|2~T6)l . □ 

Note that the angular momentum operator zd z on B has eigenvalues £ = 0, 1, 2, . . . with cor- 
responding normalized eigenfunctions fi (z) = (7r£!) -1 / 2 z e . On Sj it acts as zd z — zd z rather than 
zd z . 

Existence of a ground state 

The essential virtue of writing the Hamiltonian as in Lemma 12.11 is that the operator H B sends 
B N onto B N . Clearly (2ix)8ij, defined in (|2.14p . is a projector, in particular a bounded operator. 
Also H B commutes with the total angular momentum operator 

N 

C N = J2^9 Zj (2.20) 

3=1 

and we have 

Lemma 2.2 (Existence of a ground state). 

E B , as defined in (|2~TO|) is an eigenvalue of H B with (possibly non unique) eigenfunction F B . One 
may choose F B to have a definite angular momentum, say Lq: 

C N F B = L F B . 

Equivalently the infimum in (|1.3|) is attained. One may choose a minimizer ^ L with definite 
angular momentum Lq. 

Proof. Since H B commutes with wc may look for its ground state by looking at the joint 
spectrum of H B and Cjy, 

E B = inf inf \(F, H b F) , F <G B N ,CnF = LF\ . (2.21) 

The subspace of B N given by CnF = LF has finite dimension (it is spanned by the elementary 
symmetric polynomials of variables with total degree L) and thus the bottom of the spectrum 
of H B in this subspace is an eigenvalue. Obviously we have 

jv 

H B > N (lj + 2k) + ((w + 3fc) Lj + kLf) 

3 = 1 

and it is easy to deduce (see Lemma H . 1 1 below) that 

inf {(F, H b F) , F G B N ,C N F = LF) -> oo when L — > oo. 

The infimum over L in (|2.2ip is thus attained, say at Lq, and the bottom of a B N H B is an eigenvalue, 
with eigenfunction (a priori not unique) F B , which may be chosen with momentum Lq ■ □ 
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The yrast curve and fully correlated states 

In the case k = 0, the Hamiltonian H is simply 

H B = Nlu + luC n + gl N (2.22) 
with C at the total angular momentum (|2.20[) and In the total interaction operator 

Ijv = ^2 Sij. (2.23) 

l<i<j<JV 

This Hamiltonian has been studied extensively in the literature see e.g. the reviews [Col IViej and 
[LeSel iLSYl IRSY) for references. Its essential property is that the operators Cn and In commute. 
The lower boundary of (the convex hull of) their joint spectrum in a plot with angular momentum 
as the horizontal axis is called the yrast curve (see jRJll IRJ21 IRCJJ1 IVHRj and jLeSc , for plots 
showing its qualitative features.). As a function of the eigenvalues L of Cn the Yrast curve I{L) is 
monotonously decreasing, starting at 1(0) = CN(N — 1) and hitting zero at L = N(N — 1). The 
monotonicity follows from the observation that if a simultaneous eigenfunction of Cn and In is 
multiplied by the center of mass, (z\ + ■ ■ • + zn)/N, the interaction is unchanged while the angular 
momentum increases by one unit. 

For a given ratio ui/g the ground state of (|2.22[) (in general not unique) is determined by the 
point (s) on the yrast curve where a supporting line has slope —ui/g. For L < N the ground state of 
(|2.22p is explicitly known |PB1 ISW] while for large N and L <C N 2 a Gross-Pitaevskii description 
with an uncorrelated ground state is asymptotically correct [LSYj . For L = N(N — 1) the unique 
ground state of In with eigenvalue is the bosonic Laughlin state (|1.8|) whose wave function in 
Bn is the symmetric polynomial 

^Lau^l, • • • ,Zn) = CLau (z { - Zj) 2 (2.24) 

l<i<j"<W 

with a normalization constant CLau- More generally we have 

Lemma 2.3 (Null space of the interaction operator). 

The null space of the interaction operator (|2.14p is given by 

Ker(I N ) = {F L&U F, FeB N }. (2.25) 

Proof. Null states of In must vanish on the diagonals = Zj, hence contain the factor Y[i<j( z i~~ z j) 
by analyticity. Null state are thus of the form 

F(z 1 , . . .,z N ) = (zi - Zj)G(zi, . . .,z N ) 

l<i<j<N 

and the bosonic symmetry of the wave function F imposes that G be anti-symmetric (fermionic) 
with respect to particles exchanges 

G(. . . ,z h . . . ,z j: . . .) = -G(. ..,Zj,...,Zi,.. .), for any i,j 

Then G also has to vanish on the diagonals Zj = Zj , which leads to the form 

F(z u ...,z N )= H (z l -z J ) 2 G(z ll ... 1 z N ) (2.26) 

l<i<j<iV 

by analyticity again and proves (|2.25p . □ 

We shall call states of the form (|2.26[) fully correlated states because adding more correlations 
to them cannot decrease further the interaction energy. Their angular momentum spectrum is 
contained in L > N(N — 1) and amongst them only the Laughlin state has momentum exactly 
N(N-1). 
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Spectral gaps of the interaction operator 

Clearly, a proof that a ground state of H B almost fully lives in Kcr(ZAr) will rely on the operator 
In having a spectral gap above its ground state. However this is not known at present although 
it is widely believed to be true. It is however possible to restrict our attention to states satisfying 
bounds on their angular momentum. For example, states having too large an angular momentum 
will be proved to also have an unreasonably large potential energy as compared to that of the trial 
states we are going to construct. Once restricted to states with finite angular momentum, the 
interaction operator is a non zero operator on a finite dimensional space and thus it trivially has a 
gap. Important quantities for us will be the gaps that one obtains when restricting the interaction 
operator in such a manner : 

gap (L) := min (a (l N \{c N =L } ) \ {0}) (2.27) 

where Cn is given by (|2.20[) . The function L H> gap(L) is decreasing for the same reason that the 
yrast curve is. 

Numerical diagonalizations of the interaction operator (see e.g. |RCJJ[ iRJll IRJ21 IVHRj ) for 
small numbers of particles suggest that the unrestricted operator truly has a gap, and that it is 
attained at angular momentum N{N — 1) — 1 : 

gap (L) = g&p(N(N - 1) - 1) for any L > N(N - 1). 

It is also believed (see e.g. discussions in jLeSej ) that this gap stays of order 1 when N —> +oo. 
Proofs of these conjectures are unknown to us and hence our criteria for the ground state of p. 711 
being strongly correlated will depend on quantities whose TV-dependence is not known. If true, the 
conjectures would imply a uniform lower bound 

gap(L) > gap(iV(iV - 1) - 1) > C for any LeN 

which would simplify the conditions on the parameters to have strong correlations in the ground 
state that we give below and make explicit their dependence on N. 

2.2 Main results 

With the notation above one can easily see [LeSel Section 2.2] that when k = and 

g gap(N(N - 1) - 1) > u 

the ground state of H B is exactly given by the Laughlin stat(Q. Decreasing u any further does 
not increase the total angular momentum since the Laughlin state has the minimum angular 
momentum amongst fully correlated states. 

When k 7^ 0, the single particle part of H B is no longer a multiple of Cn as in (|2.22p . In 
particular it no longer commutes with the interaction operator and the Laughlin state is no longer 
an exact eigenstate of the full Hamiltonian. One may only hope that F B becomes almost fully 
correlated in an appropriate limit. Identifying such a limit is the main goal of the companion paper 
RSY] . Note also that the spectrum of the single particle Hamiltonian in (|2.13|) is given (up to an 
additive constant) by 

(lu + 3k)£ + k£ 2 , £eN 

with normalized eigenfunctions fi{z) = (nd)^ 1 / 2 z e . When u> < 0, which is allowed if k ^ 0, 
the single particle part thus favors states with non zero angular momentum £, in contrast with 



4 Note that gap(A r (A r — 1) — 1) > because all states with zero interaction arc of the form 1 12. 261 ) 
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the situation when u> > 0. One may thus imagine to obtain fully correlated states with angular 
momentum larger than N(N — 1). 

We now state precisely our results, starting by providing estimates on the angular momentum 
of the ground state of (|1.7|) . In particular, this will provide the reduction to states having finite 
angular momenta we were alluding to before. Since our Hamiltonian commutes with the total 
angular momentum £^r we may choose \1/ L in a definite total angular momentum sector. We 
denote L the total momentum of such a ground state : 

N 

J2( z 3 d zj ~ %<^)* L = ^o* L (2.28) 

3=1 

and we have 

Theorem 2.1 (Angular momentum estimates). 

In the limit N — > oo, w, k — > the angular momentum Lq of a ground state of H L satisfies 

1. Ifuj> -2kN 

L < 2N 2 (2.29) 

2. Ifu>< -2kN and \uj\/k < N 2 

\L -L qh \<V3N 2 , (2.30) 

where 

L qh = -^ + 0(l). (2.31) 
In particular L /L qh ->• 1 if N < \u)\/k < N2. 

3. Ifu < -2kN and \u\/k > N 2 

\L -L qh \ < VSL^N. (2.32) 

In particular Lo/L q h — > 1 i/ |w|/fc 3> ^V 2 . 

The situation described in Item 1 of the above is compatible with the Laughlin state staying the 
ground state of the Hamiltonian. On the other hand, Items 2 and 3 show that the Laughlin state is 
not the true ground state for sufficiently negative values of the ratio tu/kN: in this case the ground 
state has a momentum given to leading order by L q h (qh stands for quasi- hole again). A state 
with much larger angular momentum than Laughlin's is thus favored. As we will discuss below, 
we prove that a wave function containing a vortex at the origin in addition to the correlations of 
the Laughlin state has a lower energy than the pure Laughlin state. 

Our next theorem is a criterion for the ground state of (|2.13|) to be asymptotically fully corre- 
lated, in the sense that its projection Pkcv(x n ) ± ('^ L ) 011 the orthogonal of Ker(Zjv) (as defined in 
(jl.lOp ) vanishes in a certain limit. As anticipated above, our criteria depend on spectral gaps of 
the interaction operator. Recalling (|2.27p . let us define 

Ai = gap(2A 2 ) 

A 3 = gap(L qh + V3N 2 ) 

A 4 = gap(L qh + VSL^N), (2.33) 

where i q h is defined as in (|2.3ip . The indices in the notation correspond to the different cases 
in the following theorem. A reader willing to take for granted the conjectures about the spectral 
gaps of Xjv we discussed above may replace these quantities by fixed numbers in the following 
statements. 
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Theorem 2.2 (Criteria for strong correlations in the ground state). 

Let ^E ,L be a minimizer of (|1.3[) . We have 

||^Kcr(i JV )^* L H0 (2.34) 

in the limit N — > oo, cu,k — > if one of the following conditions holds : 
Case 1. oj > and 

(gAi)' 1 (wiV 2 + kN 3 ) -> 0. 

Case 2. -2fciV < w < and 

(5 Ai)" 1 (iV w 2 /fc + a;A^ 2 + fcA^3) _> q. 

Case 3. < ~2kN and 

(g A 3 ) _1 fciV 3 -> 0. 

Case 4. < ~2kN and 

(g A 4 ) _1 |w|JV -> 0. 

Note that, given some u; and fc depending on iV in a definite manner, one can always choose g 
so large that one of the criteria in Theorem 12.21 is satisfied in the limit N — > 00. 

We now state some energy estimates. In the cases described in Theorem 12.21 we are able to 

determine the order of magnitude of the leading order of the energy, though with unmatching 
constants. 

Theorem 2.3 (Energy bounds). 

The ground state energy _E L satisfies the following bounds: 
Cases 1 and 2 

(ujN 2 + kN 3 ) (1 - o(l)) < E L < (ujN 2 + ^kN 3 ^j (1 + o(l)). (2.35) 

Case 3 

-^(1 - (1)) < & < (-^ + '-kN 3 ^ (1 + o(l)). (2.36) 

Case 4 

-^(1 - (1)) <E L < (-^- + 1\ U \N) (1 + o(l)). (2.37) 

The proofs of these results are given in |RS Yj . The rest of the paper is devoted to their 
interpretation, using the plasma analogy that we expose in details in Section [3.11 

2.3 Discussion 

The reason for the occurrence of different cases in Theorems 12 . 21 and 12 .31 can be interpreted in view 
of the plasma analogy. Indeed, the minimization of (|1.16l) is a simple electrostatics problem and 
accurate approximations to p MF may be computed, leading to the following picture: 

• In case 1 the effective potential is increasing and it is favorable to use the pure Laughlin state 
as a trial state. We prove that in this case the mean-field density is approximately constant 
in a disc around the origin. 
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• In cases 2 to 4, the effective potential has a local maximum at the origin and a minimum 
along some circle of radius r op t(w, k). In case 2 the potential well along r = r opt (ui, k) is not 
deep enough to make it favorable to deplete the density of the trial state at the origin, and 
the Laughlin state is still preferred. 

• For the cases 3 and 4, we notice that the mean-field density g MF is well approximated by a 
profile with a maximum along a circle of radius r(rn,N). Equating r op t(w,fc) and r(m,N) 
in order that the maximum of the density coincides with the minimum of the potential we 
find an optimal choice for the phase circulation m of the giant vortex at the origin 



The Laughlin state is thus favored for to > —2kN whereas there is a tendency towards adding 
a vortex at the origin in the opposite regime. 

• The character of the density of the optimal trial state changes from an almost constant profile 
when m <C N 2 to a Gaussian profile when m ^> N 2 , corresponding to the change in the order 
of magnitude of the subleading contribution to the energy upper bound in Theorem 12. 31 that 
distinguishes case 3 from case 4. 

The fact that the mean-field density profile changes in the regime m oc N 2 is the reason why 
the estimates of Theorem 11.11 are stated differently when m <C N 2 and m ^> N 2 . The change 
in the physics is reflected by the need of a different approach to the mean-field limit in the two 
regimes. 

We thus see that the plasma analogy provides a rationale for the occurrence of different cases 
in the minimization of the energy functional describing rotating bosons in the lowest Landau level 
when the trapping potential is of the form (|1.5|) (more complicated expressions could be considered) . 
As we prove below (Section 3|) , it also allows to improve some of our energy upper bounds. The 
rest of the paper is devoted to the proof of Theorem 11.11 and to the study of the mean-field energy 
functional (|1.16j) . 

3 Quantum Hall states and the plasma analogy 

In this core section of the paper we use the interpretation of the modulus squared of fully correlated 
trial states as the Gibbs measure of a 2D Coulomb gas (one-component plasma) to compute the 
single-particle density of the trial states in the limit of large particle number. This interpretation 
has been instrumental since the first introduction of the Laughlin state [Laui ILau2] in the context 
of the fractional quantum Hall effect, see jGirj for a review. Ideas derived from the plasma analogy 
are commonly used in the literature (see e.g. |BCR1 IGRG1 ILFS1 Uanl IJLSj ). and we shall provide 
a rigorous justification to some of them. 

It is convenient to scale distances by a factor ■ 




(2.38) 



m {Z) := N N VttfNZ) 



(3.1) 



With such a rcscaling, we can recognize the Gibbs measure of a 2D Coulomb gas with temperature 
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with 



T = N 1 and a mean-field scaling in the interactions (Z^ is a normalization factor) 

( N 

p N (Z) = Zjfexp ^(-iV|z,| 2 + 2TOlog|z,|) -4^]log|^-^ 
\J= 1 «J 

= exp (-1 (|,/ - 2^ log N ) £ log I 

\ \J'=1 *<0 

= Z^-Bxp(~H N Y (3.2) 

where the Coulomb Hamiltonian Hn is defined as 

N 2 

H N (Z) := £ ^m(%) " ]v E l0g l z * " Z A ( 3 ' 3 ) 

iy rn (z) = |z| 2 -2^1og|z|. (3.4) 

This model describes N classical 2D particles located at points z%, . . . , Zn in the complex plane, 
interacting via 2D Coulomb forces and feeling the electric potential generated by a constant back- 
ground of opposite charge (the \zj\ 2 terms). When m ^ 0, the term — 2^ log \zj\ describes the 
effect of a particle of charge 2^ fixed at the origin. Our classical one-component plasma is thus 
more precisely a jellium with an additional point charge pinned at the origin. 

Note the 1/N factor in front of the interaction term: the interest of scaling the distances is 
to put us in a mean-field regime. Common wisdom about the thermodynamic limit for classical 
particles then suggests that we shall be able to extract information about pn from a limit N — > oo. 
More precisely, one should expect that pn factorizes 

p N « p® N when N -> oo, (3.5) 

for some well-chosen probability measure p <E P(K 2 ) (see Section HOI below) . in the sense that 

tiff a p® k when N -> oo and k is fixed. (3.6) 

Here we denote by pS k > the k particle density of a symmetric measure p over a Cartesian product, 
defined by integrating p^ over N — k variables: 

/i (fe) (zi,...,z fc ) := / p(z l7 . . . 7 z N )dz k+1 . . .dz N . (3.7) 



Results in this direction are given in |CLMP[ IMSl IKiel IKS] for related models. A large deviation 
result is presented in |IS/,I. Adapting these methods we could prove that p^ N ' — 1 p® k weakly as 
measures for any fixed k. For our purpose, however, quantitative estimates are needed and we thus 
use a different method. 

Our new approach to the mean-field limit of the Coulomb gas works in any scaling of the spatial 
variables but it is important to note that once a scaling has been chosen, the dependence of the 
temperature on TV is fixed. In the most convenient mean-field scaling that we have chosen above, 
the temperature is 

T = N- 1 (3.8) 
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and thus T — > when A — > oo. One could thus expect the plasma, to which we compare the 
density of our QH states, to be close to being in its ground state. These heuristic considerations 
have to be taken with caution however since the potential W m depends on both m and N. As it 
turns out there is a transition in the physics of our trial state, the dividing line being given by 
m oc TV 2 : For in <C A 2 , electrostatic effects dominate, i.e. the plasma is close to its ground state 
at zero temperature, whereas for m 3> A 2 entropy considerations dominate the physics, i.e. the 
temperature plays an important role. 

Although we shall not use this analogy, it is worth recalling the strong connection between 
Coulomb gases and Gaussian random matrices, noted first by Wigncr, see [For] for references. 
Interestingly, in the random matrix context, one also has to analyze a Coulomb gas with mean- 
field interactions and temperature of order A -1 , see [KSj . The study of the 2D Coulomb gas 
problem, or some of its generalizations, is also related to vortex systems in classical and quantum 
fluids, see e.g. |CLMPl ICYl ISS] for discussions. 



3.1 The mean-field one component plasma 

Let us now go into more details about our approach to the mean- field limit. As is well-known 
[CLMP1 IKiej and easy to prove, pn minimizes the free-energy functional (recall the temperature 
is A" 1 ) 

F N \p\ := / H N (Z)p(Z)dZ + T [ p(Z)\ogp(Z)dZ (3.9) 

JK2JV J R 2N 

amongst symmetric probability measures p G P s (IR 2Ar ). We denote by 

F N = F n [pn] = -TlogZ N (3.10) 

the minimum free energy. 

We will prove estimates relating the minimization of (|3.9[) to that of the mean-field free energy 
functional 

£ MF [p} = [ W m p + 2D(p lP ) + T f plogp (3.11) 

JR 2 JR 2 

amongst probability measures p € P(M 2 ). We denote by g MF and E MF respectively the ground 
state and the ground state (free) energy of the mean-field free energy functional. The notation 

D[p,p) = - ff p{x)\og\x-y\p{y)dxdy (3.12) 

J JR 2 xR 2 

stands for the 2D Coulomb energy. 

As usual, £ MF is obtained by restricting Fn to trial states of the form p® N , which should 
be a reasonable approximation when A is large. The main goal of this section is to justify this 
approximation by proving that 

I'M W , PN ~ Q 

in a sense to be made precise below, and with quantitative estimates. 

The transition in the physics of our trial states between dominantly electrostatic and thermal 
behaviors can be taken into account by introducing two simplified functionals that will be used to 
approximate (|3.1ip in the two different regimes. In the electrostatic regime we drop the entropy 
term and define 

£ cl [p] = f W m p + 2D{p,p) (3.13) 
Jr 2 
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with ground state g and ground state energy E whereas in the thermal regime we drop the 
electrostatic term to obtain 

£ th [p]-/ W m p + T [ plogp (3.14) 

whith ground state g th and ground state energy E th . For our computational purpose it is much 
more convenient to estimate the difference between /jff and g cl or g th because, at least with the 
relatively simple potentials W m we consider, the two latter functions are explicit, see Proposition ^. II 
below. 

The main output of the plasma analysis is the following theorem. We denote 




the minimum point of the potential W m . 

Theorem 3.1 (Plasma analogy for quantum Hall phases). 

There exists a constant C > such that we have 

1. (Mean-field limit in the electrostatic regime). For m < iV 2 and any V such that W € 
L°°(M 2 )ni 2 (ffi 2 ) 

< CN- 1 ' 2 (log JV) 1/a || W|| ia(R3) + CW" 1 / 2 !! WlU-pp). (3.15) 
Also, for some constants c, C > and for N large enough 

l$(z) < Cexp (-cN ({\z\ - r opt ) 2 - logiv)) when \\z\-r opt \ > C maxfJV 1 /^- 1 / 2 , AT 1 / 2 ). 

(3.16) 

2. (Mean-field limit in the thermal regime). For m 3> -/V 2 and any V G L°°(M 2 ) 

<C7V 1/2 m- 1/4 ||F|| L oo (R2) . (3.17) 

Moreover there exists c, C > such that 

^ ) (z)<cxp(~cN(\z\~r opt ) 2 ) when ||z|-r op t| > CN^m,- 1 '*. (3.18) 

The proof of this theorem relies on upper and lower bounds to the free energy proving that 
F N - NE MF with controlled error. Our approach can also give information on the reduced 
densities for k fixed in the limit N — s- oo see Remark 13.31 below. We state only (|3.15[) and 
(I3.17[) explicitly because they will be our main tools for estimating the energy of our quantum 
Hall trial states. What these equations say is that we can replace p^ by g cl or g th , making a 
controlled error. Note that to put them to good use in the proofs of our main results we will have 
to truncate the physical potential T4;,fc so that the norms appearing in the right-hand sides be 
finite. To estimate the error this induces we need to know that p^ has a suitable decay, which 
is the purpose that (|3.16p and (|3.18[) serve. As we will prove below, g MF decays rather fast in the 
regions where g cl and g th are small so that one can hope not to make a large error when truncating 
the physical potential. 

The proof of this result goes as follows : In Section 13.21 we first study the mean- field functional 
and prove that g MF can be approximated by g el (respectively g th ) when m <C -/V 2 (respectively 
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when m 3> iV 2 ). We also study the decay of g MF , which will provide the desired decay of /j$ 
in the electrostatic regime according to (|3.16[) . In Section 13.31 we study the mean-field limit and 
thereby relate fjr^ to g MF , which will complete the proof of ()3.15[) . Most of the arguments in this 
part apply to much more general situations than that we are directly interested in. The interested 
reader should have no difficulty in adapting our proofs to different potentials than our specific W m , 
and to other temperature regimes than T = iV" 1 . Finally Section I3T41 contains the proof of (|3.I7[) . 
As explained below we have to follow different strategies for the electrostatic and thermal regimes, 
which accounts for the different norms of V appearing in the right-hand sides of (|3.15p and ()3.17[) . 
and the different forms of the decay estimates (|3.16[) and (|3.18D . 

3.2 The mean-field functionals 

We now state several facts about the mean-field problems, some being well-known from potential 
theory (see |STj for references). 

Proposition 3.1 (The mean-field functionals). 

The following properties hold 

1. (Existence, Uniqueness). The functionals (|3.11[) , (|3.13l) and (|3.14[) each admit a unique 
minimizer among probability measures, respectively denoted g MF , g cl and g . Moreover 

< g MF < ^- a.e. inM 2 . (3.19) 

2. (Electrostatic regime). We have the explicit expression 

Q° l = 2^(o,V2) l f m = ( 3 - 2 °) 
g° l = ^~t AN ifm>0 (3.21) 



where An is the annulus of inner radius R m = ^Jm/N and outer radius i?+ = y/2 + m/N 
centered at the origin. Moreover 

D (g cl - g MF , g cl - g MF ) < CN' 1 . (3.22) 

3. (Thermal regime). We have the explicit expression 

^V) = ^cxp(-T^ m (r)) = ^z'»cxpHz| 2 ) (3.23) 
Z tn 7rm! 

where Z th is a normalization constant satisfying E = —TlogZ . Moreover, 

D (g th - g MF , g th - g MF ) < Cm- 1 ' 2 (3.24) 

for any m N' 2 and 

lk th -P MF || TV <^ (3-25) 
where ||/i|| TV = J + J stands for the total variation norm of a measure \i. 



Remark 3.1 
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1. Norms. Note that D(., .) is the square of a norm on the space of measures with total mass 
as we will see below. It is actually the square of the norm (L 2 norm in Fourier space 
with weight |fc|~ 2 ). Comparing (|3.22p and (|3.24[) one can see that g MF is better approximated 
by g° l (respectively g th ) when m <C iV 2 (respectively m ^> N 2 ). When m ^> N 2 , it becomes 
possible to use the total variation norm to estimate the difference between g MF and g th p.25[) . 
which is more convenient for practical purposes. Note that in our case g MF and g th are L 1 
functions, so their total variation norm coincides with their L 1 norm. 

2. Comparing the profiles. When m 3> N it is safe to approximate W m with its second variation 
around r opt = m 1 / 2 N~ 1 / 2 , noting that W^(r pt) = 0(1), which means that g th is roughly 
speaking a Gaussian profile centered on the minimum of W m . Obviously this is a very 
different shape from the electrostatic profile, which is constant in an annulus close to r opt . 
More important are the scales involved in the two profiles: the electrostatic density has 
its maximum of order 1 and consequently its support has a thickness of order mT l / 2 N 1 / 2 , 
whereas the thermal profile has a maximum of order m~ 1 / 2 N and thus is spread over an 
annulus of thickness N~ x / 2 to ensure normalization. 

3. Heuristics for the electrostatic/thermal transition. A good criterion for the transition, that 
can be backed with energetic considerations, is the comparison of the length scales: To favor 
the potential energy, g MF wants to be as concentrated as possible close to the minimum of 
W m , the meaning of "possible" being set by the other terms in the functional. The entropy 
and electrostatic terms are associated with different length scales, N^ 1 ^ 2 for the entropy and 
ttj-!/ 2 ^ 1 / 2 for the Coulomb term. In order to minimize the energy, the true profile g MF is 
spread over the maximum of these two length scales, i.e. on the electrostatic length scale for 
m <C ^V 2 and on the thermal length scale for m 3> ^V 2 . 

□ 

We start the proof of Theorem 13.11 by recalling some well-known lemmas that we shall use 
several times in the sequel : 

Lemma 3.1 (Positivity of relative entropies and CKP inequality). 

Let [i and v be two probability measures with /i absolutely continuous with respect to v . Then 

J Mlog^ > 0. (3.26) 

More precisely one has the Csiszdr-Kullback-Pinsker ( CKP) inequality 

l^log^ill/x-HlTv (3-27) 
Proof. A simple application of Jensen's inequality : 

since /i and v are probability measures and x i— > x log x is convex. A proof of the CKP inequality 

and some generalizations may be found in |BVj . □ 

Lemma 3.2 (Positivity properties of the 2D Coulomb energy). 

Let fi be a Radon measure over M 2 whose positive and negative parts pi + and /j,- satisfy 
\D{n+,fJ,+)\ < oo, < oo. // 

/ M = 

JR 2 
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then 

D((i,n)>0 (3.28) 

with equality if and only if fi — 0. 

Consequently, the functional \x n- D(fi, y) is strictly convex on the space P(R 2 ) of probability 
measures on R 2 . 

Proof. This is a consequence of the formula 

£W) = t^ / (f 7-^—M z )) dt 



that holds whenever J R2 ^ = 0, see [ST1 Chapter I, Lemma 1.8]. 

To see that (|3.28|) implies the claimed convexity property, pick Hi,fi2 € P(R 2 ) and notice that 

1,.1.,/llllNl, 

-D{in, m) + -D(ji 2 ,^2) — D I -/ii + -^ 2 , -/ii + -/i 2 = 7 £> (Mi - M2, Mi - M2) > 



since J" R2 /Mi = J* R2 // 2 = 1- □ 

Lemma 3.3 (Newton's theorem). 

For a measure \x let 

K{x) = - J log\x - y\n(dy) (3.29) 
be the potential generated by \i. If fi is radial then 

h^x) = -\og\x\ ( ii(dy)- [ log\y\fx(dy). (3.30) 
J\v\<\*\ -%I>M 

Proof. Simply reproduce the proof of the corresponding result in 3D, see |LL[ Theorem 9.7]. □ 



Proof of Theorem \3.1\ 

Step 1. (Existence, Uniqueness) . The existence part is standard material, as are the following 
Euler-Lagrange equations 

ttjHF + W m + N- 1 log(g MF ) = E MF + 2D(g MF , g MF ) on R 2 (3.31) 
W m + N- 1 log((? th ) = E th on R 2 . (3.32) 

Strictly speaking, studying the variations of £ MF and £ th only tells that the above equations hold 
on the support of g MF , respectively g th . However, as we prove below, g MF and g th are strictly 
positive a.c. which is why the above equations hold on R 2 . 

On the contrary, g cl has compact support (see below) , which makes the Euler-Lagrange equation 
a little bit more subtle ISTI : 



4/i e oi + W m = E cl + 2D(g c \ g cl ) on supp( £ - cl ) (3.33) 
4/i e d + W m > E cl + 2D{g c \ g cl ) on R 2 \ supp(# el ). (3.34) 

Here h g MF and h„ei are defined as in (|3.29[) . The value of the constants on the right-hand sides 
of the Euler-Lagrange equations is evaluated by multiplying the equations by g MF , g th and g el 
respectively, and integrating. The formula (|3.23[) is a direct consequence of p.32[) once one knows 
that g th > a.e. and thus that (|3.32[) holds on the whole space. 
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The fact that we have equality on the whole of M 2 in (|3.3ip and (|3.32[) follows from the fact 
that supp(g MF ) = supp(g th ) = M 2 , which is probably a point that deserves a little discussion. We 
follow an argument from [Ner| (proof of Proposition 15 therein) : Suppose for contradiction that 
supp(g MF ) c contains a set S of nonzero Lebesgue measure. Consider the trial state 

g MF + et s 

P : 



l + e\S\ 

for some e small enough. Evaluating £ cl [p] is easy and we find that for small enough e there is a 
constant C > such that (we consider TV, m and T as fixed here) 

£ cl [p] <£ cl [g MF ]+Ce. 

To compute the entropy of p, the key point to notice is that since g MF and els have disjoint 
supports we have 

(g MF + el s ) log (g MF + el s ) = [ g MF log g MF + e [ l s log (el s ) 

Jr 2 Jul 2 

< [ g MF log g MF + Ce(l + logs), 
Jr 2 

and thus 

£ MF [p] < £ MF [g MF ] + Ce (1 + logs) < E MF 

for e small enough, which is a contradiction. The proof that also supp(p th ) = K 2 is identical. 
To see that g MF < (2tt)~ 1 a.e. in M 2 , we take the Laplacian of (|3.31[) and obtain 

Ag MF |Vg MF | 2 o mf ,,4^ 

- T ~pM^ +T ( g MF) 2 + 8?r g -4+-^-5 -0 

which implies 

-TAg MF + 8ir(g MF ) 2 - 4g MF < 0. 
At any local maximum of g MF we have Ap MF < and thus 

which proves the claim. 

The uniqueness of the minimizer of £ MF follows from the strict convexity of the functional. 
Quantitatively, we have the following stability identity: For any probability measure p that we 
write as p — g MF + v, 



£ MF [p}= [ W m g MF +f W m v + 2D(g MF ,g MF ) + 4D( t 

JR 2 JR 2 



g MF ,v)+2D(v, v) 



Tl (g MF + v)\og(g MF + v) 



= 1 W m g MF + 2D(g MF ,g MF ) + T f (g MF ) log (g MF ) + 2D{u,v) 

JR 2 JR 2 

+ T [ ((g MF + v)\og(g MF + v) -v\ogg MF - g MF \ogg MF ) 
Jr 2 

= E MF + 2D{v, v)+T f (g MF + v) log (^^) (3.35) 

Jr 2 V Q J 
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where we have used (|3.3ip and the fact that J R2 v = to go to the second line. This yields 
uniqueness for g MF using Lemmas 13.11 and 13.21 Note for later use that this also proves stability of 
the minimizer in H^ 1 and TV norms. 

Step 2. (The electrostatic profile). The proof of ([3~zT))) and ([3~zTj) relies on ([535]l . We start 
with the easiest case m — where we obtain the circle law for the Ginibre ensemble (see |KSj 
for example). A proof of (|3.20[) may be found in |STj . we give details for the convenience of the 
reader. Taking the Laplacian of (|3.33[) we have 

g cl = — on supp( £ - cl ). (3.36) 
2tt 

By uniqueness of the minimizer of £ el and the radiality of W m , g o1 must be radial. Its support is 
thus in any case a union of some annuli (counting a disc as an annulus with inner radius zero). 
Suppose for contradiction that the support is not a disc. Then there is a nonempty annulus A (or 
a disc centered at the origin as a special case) in the complement of the support and enclosed by 
the support. The potential 

tp = 4/ V , + W m - E cl - 2D(g cl , g cl ) (3.37) 

satisfies tp = on dA by (|3.33[) and — Atp = — 2 < in A since g cl = there by definition. We 
deduce by the maximum principle that <p < and thus tp = in A because by (|3.34[) we already 
know that tp > 0. Taking the Laplacian of the equation tp = on A we would conclude as before 
that g cl = (2tt)~ 1 in A which is a contradiction with the fact that ^4(~lsupp(p el ) = 0. We conclude 
that the support of g el must be a disc centered at the origin, which implies f|3 . 20[) via (|3 . 36[) and 
the normalization constraint. 

The result in the case m > follows from the same kind of arguments. We have from (|3.33l) 

el 1 m /oh 

Z7T ZiV 

and since g cl > 0, we deduce that ^ supp(g cl ). By the same maximum principle argument as 
above we deduce that supp(g> cl ) is an annulus and that (|3.36p also holds in the case m > 0. We 
consider again tp as defined in (|3.37|) . By (|3.33l) . it must be that tp = on supp(g el ). In particular 
it is constant there. On the other hand, 

p = -E cl - 2D(g c \ g cl ) + tp = const + tp 

is a constant plus the potential tp — Ah g d + W m generated by 

• a constant background of charge density — (2tt)~ 1 (coming from the \z\ 2 term in W m ) 

• a point charge of strength ^ located at the origin (coming from the 2^ log \z\ term in W m ) 

• the charge density 4g> o1 , equal to (2ir)~ 1 in an annulus of radii say i?~ and i?+ and elsewhere. 

Using Newton's theorem, Lemma 13.31 this potential, evaluated at any i?~ < r < i?+ equals that 
generated by a point charge of strength =^ — 2i?~ 2 located at the origin plus another constant 
(the constant background and the charge density Ag cl cancel each other in supp(p o1 )). The only 
possibility for tp, and therefore tp, to be constant in supp(p cl ) is then to have 
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so that the effective potential generated on supp(p cl ) by the constant background in B(0, i?~ ) is 
canceled (screened) by the point charge sitting at the origin. Using the normalization of g el and 
(I3.36[) . we compute i?+ and (|3.21[) follows. 

Step 3. (Electrostatic regime). We now turn to the proof of (|3.22[) . First, taking g el as a trial 
state for £ MF we have 

E MF < £ol + T f g ei log ^cl 



which, in view of (|3.20p and (|3.21[) . yields 

e mf < ^ci _ log ( 27r )T. (3.38) 

Next, we denote 

p (z) = 7r _1 exp(-|z| 2 ) 
and note that J R2 po = 1 . Using Lemma 13.11 we have 

e mf = 2D(g MF ,g MF )+ f W mQ MF +T f g MF \ogp + Tf g MF log ^- 

Jr 2 Jm 2 Jr 2 Po 

> 2D(g MF ,g MF )+ f (W m T\z\ 2 ) g MF - log(n)T. (3.39) 

Jr 2 



Now, the functional 



£ C \ P ] :=2D(p,p)+ [ (W m -T\z\ 2 )p 

JR 2 



is of the same type as £ el , at least if N is large enough for the —T\z\ 2 term in the above to be smaller 
than the \z\ 2 term in W m (recall that we have T = N^ 1 ). We denote E cl and g cl respectively the 
ground-state energy and the minimizer of £ . Using the Euler-Lagrange equation satisfied by g cl , 
similar to p.33p . and a computation analogous to (|3.35[) wc obtain for any measure p 

£ c \p] > E cl + 2D(p- g c \ p - g cl ) . (3.40) 

Note that in this case there is no entropy term in the analogue of (|3.35[) and the second line becomes 
an inequality because g el has compact support and thus it satisfies an Euler-Lagrange equation 
only on its support. Outside of the support we have an inequality as in (|3.34p . We deduce from 
(pHp and [ESQ) that 

E MF > E el + 2D (~cl _ ? MF ; -el _ g MF) _ CT ( g M) 

Using the explicit expression for g e , similar to (|3.21[) . it is then not difficult to see that also 

E cl = E cl + ( T ) 

and thus p.38p and (|3.41[) combine to give 

D (g MF - g°\ g MF - g cl ) < CT 
and we deduce that (|3.22p holds by noting that also 

D (g cl - g c \ g cl - g cl ) < CT 
which can be proved easily, inspecting the explicit expressions for g el and g . 



2G 

Step 4- (Thermal regime). Taking g th as a trial state for £ th we have 

e mf < E th + 2D(g th , g th ). (3.42) 

To obtain a lower bound we write 

D(g MF , g MF ) = D(g MF - g th , g MF - g th ) - D(g th , g th ) + 2D(g th , g MF ) 

and deduce 

77.MF ^ cthr MFi 1 orV MF th MF th\ on/' th th\ /o ao\ 

E > t- [g \+2L){g —g ,g —q ) - 2V(g ,g ) (3.43) 
where (we denote h„th the potential associated to g th ) 



£ th [p] := f (W m + 4fc B *) p + T f plogp 
Jm 2 Jr 2 



with ground state g th and ground state energy E th . Of course 

g th = ± exp {-T-\W m + Ah gt ,)) (3.44) 

for some normalization constant Z th satisfying 

E th = -TlogZ th . (3.45) 



Now, h g ti> is radial and satisfies on R 2 



-A/i e th = g th > 0. 



Integrating this equation over S(0, r) and using Stokes' theorem we deduce 



ih 



2TTrd r h„th(r) = — J g 

JB(0,r) 

and thus |V/i e t h (i)| < CY" 1 for any t £ B (r/2, 2r). From this we deduce the estimate 

\h g th(r) - fo e th(r opt )| < Cr~p t \r - r opt | for any r E B Qr opt , 2r opt ^ . (3.46) 

Since on the other hand W m — min W m grows as C(r — r opt ) 2 close to r opt , it is easy to deduce from 
()3.44j) that g th is exponentially small in the region where g th is, that is for \r — r opt | ^S> N~ 1/>2 . 
Indeed, note that in this region 

\W m - minW m | oc C\r - r opt | 2 > {N/m) 1/2 \r - r opt | oc r~ p \\r - r opt | oc \h g t h (r) - h g t b (r opt )\ 

provided m ^> N 2 . Simple estimates then show that 

E th = -TlogZ th = -Tlog J ^ exp (— T _1 Wm) exp (<ih gth (r opt ) + O^- 1 / 2 ))) 

= -Tlog [f exp (-rV m )j + 4/ leth (r opt ) + O^ti" 1 / 2 ) 

= £ th + 4h gth (ropt) + 0(m~ 1/2 ). (3.47) 
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On the other hand, similar considerations based on (|3 . 23[) and (|3.46[) lead to 

D(g th , Q th )= [ h etb g th = h^(r opt ) + 0(m^ 2 ) (3.48) 



where the last term is r~ p \ = \J N jra times the length scale of g th (we are basically saying that g th 
resembles a delta function concentrated along the circle of radius r opt ). Coming back to (|3.43|) , 
using (|3.47p and (|3.48|) , we have thus proved that for m ^> TV 2 

E MF > E th + 2D(g th , g th ) + 2D(g MF - g th , g MF - g th ) + 0{m~ 1 / 2 ), 

which we combine with (|3.42[) to obtain (|3 . 24|) . 

To prove (|3.25p we go back to (|3.43p again and note that for any p E P(M 2 ) we have 



»g^>^ h + |||p-^ th ||TV 



as a consequence of the explicit expression of g th and the CKP inequality (|3.27|) . Combining this 
with the considerations above, our lower bound can be improved to 

P MF ^ cith , on/ th th\ . on/ MF th MF th\ , T II MF ~th|| 2 , r^/-l/2\ 

h, > hi + 2D{g ,g ) + 2D(g -g ,g —g ) + — \\g — g || TV + (J(m 1 ). 

Combining with the upper bound p.42[) we deduce 

\\g MF - £ th || TV < CT^m- 1 ^ = CWV2 m -i/4 

and (|3.25p follows by estimating the difference between g th and g th , using (|3.46[) and the explicit 
expressions (|3.23[) and (|3.44[) . 

□ 

As announced, the proof of (|3.16[) requires an estimate of the decay of g MF . This is the content 
of the following 

Proposition 3.2 (Decay of the mean-field density). 

There exists a C > such that for any r £ R satisfying \r — r opt | > Cmax(A rl / 2 m -1 / 2 , N^ 1 / 2 ) 

g MF (r) < Ccxp(-C7V(r - r opt ) 2 ). (3.49) 
Proof. Step 1. We start by proving that the potential 

/Vi_ e MF = V 1 - V IF = 27r (- A ) _1 (/ - Q MF ) ( 3 - 50 ) 
is in £P(R 2 ) and that one has the bound 

|VV_ s m P | 2 < CT. (3.51) 

Let us denote / the Fourier transform of a function /. We have (formally for the moment) 

h e ,_ e ^{k) = ^{g MF {k)~g c \k)). 
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Now, since both g and g are uniformly bounded in L 1 and in L°°, they also are uniformly 
bounded in L 2 , which implies that g MF — g cl is uniformly bounded in L 2 . On the other hand, since 
J 2 W m g MF and J R2 W m g el are finite we deduce 

/ \x\g MF (x)dx < oo, / \x\g el {x)dx < oo 
Jm. 2 Jr 2 

which implies that V (g MF - g cl ) £ L°°(IR 2 ). It remains to recall that 

-MFm\ *ol/n\ / I MF cl\ n 

g (0) — g (0) = J 2 \Q ~ Q ) = 
to deduce that, for |A;| small enough 

|£ MF (fc)-£ cl (fc)| < C|fc|. 
Together with the fact that g MF - g cl e £ 2 (R 2 ) this implies that 

/ |VV- e M F | 2 =47r 2 / J_|£ MF (fc)-£ el (fc)f <oo. 

./k 2 7 R 2 |fc| 2 

We can thus justify the integration by parts leading to 

f |VV- e M F | 2 = [ -AV- r (f d -? MF )=D(/-f MF ,^ L f MF ) 

JR 2 JK 2 

and (j33T|) then follows from ([3722]) . 

S^ep 2. We now claim that for any r6l 

4Vi_ eMF (r) < S cl + 2D(g cl , g cl ) - E MF - 2D(g MF , g MF ) 

+ c ( T + Tl/2 (r-r op ^Hr + r 

\ min(r opt ,r) J 

First note that taking the difference of equations (|3.3ip and (|3.33[) we obtain 

4/ ie o 1 _ eMF (r opt ) = E ci + D(q*\ g ci ) - E MF D(g MF , g MF ) + Tlog g MF (r opt ) 
< E cl + D(g c \ g cl ) - E MF - D(g MF , g MF ) + CT 

by using (|3.19[) . Using radiality we then have 
V'-s UF ( r ) 

< 

where we use (|3.51[) and the Cauchy-Schwarz inequality. This proves (|3.52l) . 
Step 3. Using (pPTj) we have 




h g ei_ e uv{r opt ) + / h'i_MF{t)dt 



tr-o» 

opt 



(v) + ^^(j VVl - MF|<(it 

h g ei_ g uw (r opt ) 



pt 

1/2 



CT 1 / 2 r^_rl 2l 
min(r op t,r) 1 2 2 
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Inserting (|3.52[) and using h ge i_ g MF = h ge \ — h g uF we deduce 
1 



g MF (r) < exp 



T 



(E* + 2D(g c \ Q cl ) - W m {r) - ih^(r)) 



x exp 



i/ a (rapt + r) 1/2 (r qP t-r)V ii 
min(r opt ,r) 



(3.53) 



Then the exponential fall-off in Q3.49P is provided by the decay of 

E cl + 2D(g c \ g cl ) - W m (r) - Ah g «(r) (3.54) 
away from the support of g cl as a consequence of (|3.34[) . More precisely, since (|3.33p implies 

E cl + 2D{g c \ g cl ) = W m (R+) + 4h g *(R+), 
we see that (|3.54p decays as 

W m (R+) + 4/ V i(i?+) - 4VM ~ W m (r) 

for r > i?+ (we only detail this case, the proof is the same for the region r < i?~). Reasoning as 
when proving (|3.46p . one easily sees that 

\4h g *(R+) - 4VM| < Cr- p \\r - R+\ < Cm^N^r - R+\. 

On the other hand, approximating W m by its second variation around r op t we have 

W m (R+ ) - W m (r) a C(r - r opt ) 2 - C(R+ - r opt ) 2 = -C ((r - R+ ? + 2(r - i?+ )(i?+ - r opt )) 

and since — r op t| = 0(m~ 1 / 2 N 1 / 2 ) we deduce that 

£ cl + 2D(g c \ g cl ) - W m (r) - 4h g *(r) < -C(r - i?+) 2 



(3.55) 



for r > i?+ + c m 1 / 2 N 1 / 2 with well-chosen c and C. 

This decay compensates for the other terms in (|3.53[) as soon as 



\r - Rt. I » max 



T, (miV)" 1 / 6 . 



(3.56) 



Note that the first length scale in the max above is of the order of magnitude of the thickness of 
the support of g cl . The second length scale accounts for the mass spreading due to the entropy 
term and the third is associated with the "error" term 



,1/2 fropt +r) 1 / 2 (r pt-r) 1 / 2 = (r opt + r)V 2 (r opt - rf 
min(r opt , r) N 1 / 2 min(r opt , r) 



/2 



Indeed, for example when r is sufficiently close to r opt , 

(r ° P Vv2 )1/2( r° Pt "f /2 « (^o Pt )- 1/2 |r - R + J 1/2 oc (mN)-yy R+\V\ 
N L ' Z mm(r opt , r) 



using the approximations r f=a r opt and \r — r opt | fts \r — i?+ | and recalling that r opt = y/m/N. To 
see that the third term in the right-hand side of (|3.56[) is really an error, recall that we use (|3.49p 
when T = N^ 1 . There are then two cases 
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• m <C N 2 , which implies 

IN 



171 

i.e. the electrostatic length scale dominates the error and the entropic length 
• m N 2 , in which case 

\/T»(toA0~ 1/6 » \/-> 

V m 

i.e. the entropic length dominates the error and the electrostatic length. 

In both cases the error term is dominated by either the electrostatic or the entropic term. Recalling 
that |r opt - R+\ < CN 1 /' 2 m~ 1 / 2 we thus have proved that flOg) holds if 

m < and \r - r opt | > \ — 

V m 

or 

m^> N 2 and \r - r opt | > N~ 1/2 , 
which is the desired result. □ 

3.3 Thermodynamic limit in the electrostatic regime 

We now turn to the study of the large N limit of (|3.9|) . This is a rather classical question, especially 
since we are in a mean-field scaling. A line of attack for this kind of statistical mechanics problems 
has been pioneered in MS for regular intcrparticlc interactions and then carried on independently 
in [KiellCLMP] in the case of logarithmic interactions. These works deal with the regime T = 0(1) 
in which all three terms in p. lip are of the same order of magnitude. They also consider the case 
of negative temperature which is more involved but irrelevant in our context. In connection with 
several ensembles of random matrices, KS extends this approach to the regime T oc iV" 1 and 
more general Hamiltonians. 

Common to these approaches is a compactness argument, which does not lead to quantitative 
estimates on the precision of the mean-field approximation. As far as we know, it is only very 
recently |SSj that constructive estimates have been obtained, in the specific case of the 2D Coulomb 
gas that we consider here. We could employ some of these estimates (in particular Theorem 3 
therein) in our context but they would not be sufficient. We thus prefer to implement a new 
method that gives different estimates, more suited to our purpose, with a simpler proof. 

Theorem 3.2 (Mean field limit for 2D Coulomb gases). 

There exists a constant C > such that, for N large enough, we have 

1. Upper bound. 

F N < NE MF - D(g MF , g MF ) (3.57) 

2. Lower bound. 

F N > NE MF - ^ - C (3.58) 

3. Estimate on the first marginal of the Gibbs measure. 
For any V e H^M 2 ) with W <G L°°(R 2 ) we have 



1 /2 

(/$ - 9 MF ) V < C O^P) ||VF|| i2(R = ) +CiV- 1 / 2 ||VF|Uoo (E2) . (3.59) 
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Remark 3.2 

1. This result is not limited to the particular type of potential W m we consider in the plasma 
analogy. Our method can accommodate any potential as soon as the associated mean-field 
functional and its minimizers are reasonably well-behaved. Note also that the estimates apply 
to any T, the two most interesting regimes being T = 0(1) (more natural from the Coulomb 
gas point of view) and T = 0(N~ 1 ) (in relation with random matrices and quantum Hall 
phases). 

2. In [SSJ a different scaling convention is used, their j3 being given as /3 = N~ 1 T~ 1 in our 
units. The approach therein is limited to f3 > 0(1) and a fixed potential W m . It shows 
that the — log(iV)/2 term gives exactly the second order correction in the regime (3 > O(l) 
(i.e. T cx iV -1 ) with fixed potential. It moreover gives the exact third order correction to 
the free energy Fn in the limit /? — > oo, that is T -C N -1 (see |SS1 Theorem 1]). This is 
more difficult, since this connection is related to a Coulombian renormalized energy, whose 
definition is rather complex. The approach we develop here allows to recover the lower bound 
on i*jv in the regime where (3 is bounded above, i.e. T cx iV -1 . To see this, compare Theorem 
13.21 with Theorem 1 in |SSj . keeping in mind that when T cx .ZV -1 , the entropy term in E MF 
is a lower order correction : E MF = E cl + 0(N~ 1 ). 

3. In the regime T oc N~ x , we could most likely adapt arguments from jSSi Section 4 and 7] in 
order to construct a trial state capturing exactly the — log(iV)/2 in the upper bound. This 
would be rather technical because, in contrast with what is assumed in ISSj . our mean-field 
densities do depend on N when m ^= 0. We thus content ourselves with a non optimal upper 
bound that has only a marginal impact on our main theorems. 

□ 

We need two classical lemmas. The first is the 2D version of the so-called Onsager lemma (see 
e.g. |LiSel Lemma 6.1]): 

Lemma 3.4 (2D Onsager lemma). 

Let p be a radial probability measure on R 2 . Denote, for some I > 
We have, for any p such that J" R2 p = N and any (xi, . . . , ccjv) £ M. 2N 

N N \ N N 

y^i>. -i> y,i>, -D(p, P )+2j2D(p,fi Xi )-j2 D (f ji 

Xi ) ftxi )• 

i—1 i—l / i—l i—1 

(3.60) 

Proof. By Newton's theorem, Lemma 13.31 the radiality of p, implies 

with (by the way) equality if p is supported in the unit disc and minij \xi — Xj\ < I. The rhs of 
the above is equal to the rhs of (|3.60p . which proves the lemma. Note that 

/ N N \ 

D I p - ^ Mx», P - ^2 Hxi j > 

V i=l i=l / 



E 

i+3 



log 



> D 
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using Lemma 13.21 because by assumption 



i=l / 



□ 



The first term in the right-hand side of (|3.60p is usually dropped to obtain a convenient lower 
bound to the Coulomb Hamiltonian. The core of our argument consists in obtaining a bound on its 
expectation value in the Gibbs measure from our upper and lower bounds to the free energy and 
using it to control the fluctuations around the mean field density. As we shall prove, this allows 
to obtain estimates on the marginals of the Gibbs measure and in particular (|3.59|) . 

The functions p Xi in Lemma 13.41 should be thought of as unit charges smeared over small balls 
that we use to replace the point charges of the Coulomb gas. This is essential in our approach but 
has some cost that we quantify in the next lemma, which is an adaptation of a well-known lemma 
used by Lieb and Oxford (cf |LiSe[ Chapter 6] for references). 

Lemma 3.5 (The cost of smearing out charges). 

Here p denotes the normalized (in L 1 ) characteristic function of the disc of radius I. For any 
p e L°°(M 2 ) 

\D( Pl 6 Xt -n Xz )\<Cl 2 \\ P \\ Loo (3.61) 
Proof. We denote h v the potential associated to a charge distribution v. By Newton's theorem 

in M. 2 \ B(Q, I). We are thus left with computing 

p(K*i - h s* t ) 

5(0,0 



which is easily found to be equal to 



/ 1 r r 2 -l 2 \ 

The second term is easily bounded by a CI 2 ||p|| ioo while the first term is proportional to (adapting 
the analogous computation in 3D, cf [Lie] and jLiSel Chapter 6]) 

' Fir) 

-^-dr (3.62) 
Jo r 

with 

i r 27T r 

F(r) := — / / p{s,e)sdsd6. 
27T Jo Jo 

Using 

l^(r)| < ML- y, 

we can bound (|3.62[) and complete the proof of the lemma. □ 
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We can now proceed to the 
Proof of Theorem \3.2\ 

The upper bound is proved by taking the trial state ( g MF )® N in ([gTty . (f3~TU|) . The 
— Z?(g MF , g MF ) error term comes from the fact that there are N(N — 1) pairs of particles, to 
be divided by the 1/N mean-field scaling factor. 

For the lower bound we use Onsager's lemma with p = Ng MF , fi the normalized indicative 
function of the unit ball and I = N^ 1 / 2 to obtain 

N / N N \ 

H N (x 1) ...,x s )>J2W m (xi) + ^DlN(^-Y,te t ,N(?*-Y,l i ** 

i=l \ i=l i=l J 

N N 

- 2ND(g MF , g MF ) + 4]T D(g MF , - 2 £ ±D(» Xi , ^ ) (3.63) 

i=l i=l 

and we will drop the second term of the right-hand side, which is positive, for the moment. 
We now invoke Lemma [33] and use (|3.19[) to claim that 

D(g MF ,fi Xz ) = D(g MF , S Xi ) + 0(jV _1 ) = hp* fa) + C^iV" 1 ) 

which turns into 

4D(g MF ^ Xt ) = E MF + 2D(g MF , g MF ) - W m { Xl ) - T log g MF ( Xl ) + O^ 1 ) (3.64) 
thanks to the variational equation (|3.31[) . Inserting into p.63[) we obtain 

N N 

H N (x u ...,x N ) > NE UF -Tj2loge UF (xi)--J2D(v Xi ,n Xi )-C 

t=i i=i 



NE UF -T^2logQ UF ( Xi ) - D(^,f,) +log(-y) - C (3.65) 



where the second line follows from a simple computation (recall that fj, Xi is a unit charge smeared 
over the ball B(xt, N^ 1 / 2 )). There only remains to use (|3.65[) to compute a lower bound to the 
free energy of /Ujv : 

FnM > NE MF - C - i logiV + T [ mjv ( log mat - V log g MF ( Xl ) ] dan... dx N 

2 J R 2N \ ^! J 



1 i„„ at , m f .. i„„. ( 



N E - C - - logiV + T / mjv lo. 



2 b ./»2« ^ JV & V g MF ® N 



which ends the proof of the lower bound since the last term is the relative entropy of [In with 
respect to (g MF )® N and hence is positive by Lemma \'3. II 

Going back to (13.63P and retrieving the positive term we had discarded we see that our upper 
and lower bounds to the energy imply the additional estimate 

/ n N { Xl ,...,x N )D [Nq mf - Yu Xi ,Ng MF - V/Xx, ) d Xl ...dx N <CNlogN. (3.66) 

This bound quantifies how close /ijv is to (q mf )® n . In particular it implies (|3.59p . a fact that we 
state as a lemma: 
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Lemma 3.6 (Using the Onsager term). 

For any V : R 2 i— > R regular enough and any symmetric probability measure hn € ^((R 2 )^) 



(^(z)^ MF (z)) 



< cat 1 ' 2 IIVVII 



CN' 1 II Wl 



L 2 



' . / iV W \ \ !/ 2 

y R2N ^(Z)Z?(^ MF -^^ t ,^ MF -^^JJ . (3.67) 



Proof. Given a one-body potential V we start with 

Hn(xi, . . . ,x n 



) (j2 V ^ ~ N j v , V 8 MF ^j dx i--- dx N, (3. 



which follows from the symmetry and normalization of /xjv- Next, note that (with \x Xi defined as 
above) 

N f. N f. N 

(3.69) 



E%)= f V T, S *= f vY,^+0[VN\\WV\\ Lao ) 

1=1 JR2 i=l JR2 i=l 

where we use that fi Xi is a unit charge smeared over a ball of radius N~ 1 / 2 . Thus 



JR 2 V / JV J K 21V 



cfcci . . .dx 



N 



O N' 1 / 2 II Wl 



(3.70) 



Then 



R 2 

< C||W| 



L 2 (R 2 ) 



1/2 



(iV AT \ 



1/2 



(3.71) 



where 



AT 



/i = ^ e MF_ E£i ^ = 2tt (-A) 1 A^ MF - J2 



To justify these computations we argue as in Step 1 of the proof of Proposition ^. 2[ using that the 
Coulomb kernel is the Green function of —A in R 2 and that 



Nig 



MF 



N \ 



35 



On the other hand 



dX(j, N (xi, . . .,x N ) 



1/2 



< I / dX^i N (xi 1 ...,x N 
Ir 2N 



by the Cauchy-Schwarz inequality and using the normalization of fiN- Using (|3.71[) . we thus have 
dXfi N (x!, . . .,x N ) 



J^V(z) ^ MF (z)-fj^(z)^) d: 



1/2 



<C\\\7V\\ L2{R2) 



that we insert into (|3.70|) to conclude the proof. 

Our final estimate p.59[) follows by combining p.66[) and (|3.67p . 



□ 



□ 



Remark 3.3 [Estimates for higher-order marginals] 
Our approach can give estimates on not only the first marginal (iff but on any reduced density 
/4y , provided k is suitably small as compared to ./V (for example, k fixed when TV — > oo). This can 
be seen to be a consequence of our main technical estimate (|3.66[) as follows: Consider for example 
a smooth 2-body potential V 2 {z\,z 2 ) and evaluate 



/ ^n v 2 = Tt2 / ^N(xi,...,x N )V 2 (x t ,x J )dX 

= / Hn(xi,...,x n ) ( V 2 (zi, z 2 )\i Xi (zi)n Xj (z2)dzi dz 2 ) dX 

N i<7T<n J ^ 2N \J* / 

+ O (V 1 / 2 sup ||VV a (y, OIL-) + O (V 1 / 2 sup ||W 2 (., y)\\ L ^j 
as in (glal) . Then 



V] / ^(^l, z 2 )n Xt (z 1 )fi Xj {z 2 )dz 1 dz 2 - iV 2 / 

itj JR 4 JR 



)dzidz 2 - N z / Vfefo, z 2 )^ (zOe*** (z 2 )d Zl dz : 



N 



1/2 



< CN ( S up\\X7V 2 (y,.)\\ L . 2 +sup\\VV 2 (.,y)\\ L2 ) D ( N g MF - V Mxi , TV g MF - V Mxi ] 
using the fact that iV _1 X^=i Mxj and £ MF arc normalized in L 1 and arguing as in (|3.71[) . There 



3G 



only remains to use the Cauchy-Schwarz inequality as before to obtain an estimate of the form 



V2 (j$ (g MF r 2 ) 



< CN-^ilogN) 1 ' 2 sup ||W 2 (y, .)\\ La + sup || W 2 (, y)|| i2 
\ y y 

+ CN~ l > 2 (sup || W 2 (y, OILoo + sup \\VV 2 (.,y)\\ L 



which generalizes (|3.59[) . Estimates for fJ^\k > 2 follow along the same lines. 



□ 



Remark 3.4 [More information when T is larger] 

Note that we have dropped one term that we could have estimated in the proof above, namely we 
have 



fi N log ■ 



,MF 19 N 



< CT~ log N. 



(3.72) 



/R2N Q» 

This estimate is useless in the regime T = iV _1 which interests us most but can become interesting 
when T is larger, in particular in the somehow more natural case where T is fixed. 
Indeed, using subadditivity of entropy (see e.g. [Kiel Proposition 1]) 

UN log > N log (i$ 



R 2 



and the fact that 



(I3~72|) implies 



This can be turned into 



^log(^) 



MF\®N 



N 



JR 2 



i0 S mf - C TN ' 



tv - V TN 



(3.73) 



where 



TV 



stands for the total variation norm, thanks to the Csiszar-Kullback-Pinsker inequality 



recalled in Lemma \3. II As already mentioned, (|3.T3[) becomes interesting only for relatively large 
temperature (take e.g. T fixed independently of N), in which case it is a somehow better estimate 
than ([539)1 ■ 

Our approach to the mean-field limit in the thermal regime, presented in the next section, is 
based on this kind of considerations. 

□ 

We now conclude the proof of Theorem 13.11 Item 1. Our main estimate (|3.15[) follows by 
combining (|3.59[) and (|3.22[) . Indeed, using the Cauchy-Schwarz inequality in Fourier space as 
before, 

1 < cd^ -rt^m^^. 



(q mf - e el ) V 



We conclude this section by giving the proof of (|3.16[) . We use the explicit expression (|3.2[) and 
the lower bound p.65[) to obtain 



N 



MF 



J =1 



log TV C 
2T + T 
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On the other hand, we recall that 

F N = -T\ogZ Nl 

so, using (|3.57[) we have 



and wc conclude that 



Z N > cxp [ --E 



( N 1 AM N 

J2 log 9 MF (*j ) + J < II ^? MF (** ) ex P( C7V lo S ^) 

which implies (|3.16D after an integration over Z2, ■ ■ ■ , Zn and the use of (|3.49j) . 
3.4 Thermodynamic limit in the thermal regime 

Wc now turn to the proof of (|3.17|) . In contrast to what we did in the preceding section to prove 
(|3.15[) we do not work in two steps, first relating /jff to g MF and then g MF to g th . The reason is that 
the errors in the energy estimates produced by the use of Lemmas 13.41 and 13.51 are not sufficiently 
small compared to the difference between g MF and g th proved in (|3.24p . One can however rely on 
a different strategy, treating the two-body term in (|3.9I) as a perturbation of the one-body part. 
Indeed, since we were able to prove that g MF is close to g th when m 3> TV 2 , we already have an 
indication that the two-body Coulomb term is not so important in this regime. 

Proof of (I3~T71) . 

It is clearer to work in variables where the one-body potential takes its minimum when r = 1, 
with value 0. Scaling distances by a factor (m/iV) 1 / 2 it is equivalent to minimize the free-energy 
functional 

F N [lj]:= [ H N (Z)n(Z)dZ + — [ n(Z) \o gf i{Z)dZ (3.74) 

JR2JV 171 J R 2N 

with the rescaled Coulomb Hamiltonian 

N 2 

H N (Z):=J2 W m (zj ) - - J2 lo S I* - *j I ( 3J5 ) 

3=1 »#j 

and 

^ m (^) = |z| 2 -21og|z|-l. (3.76) 

Note that we have taken advantage of the normalization of /i to subtract 1 and have W m (r) > 
W m (l) = for any r € M. We can rewrite the free energy as 

77 r i Ar f f W m (x) , W m (y) 2(JV-1) , A (2) 

J r N [n\=N — 1 log \x - y /i l ' (x, y jdxdy 

J R 4 V 2 2 m J 

+ — [ n(Z)\ogn{Z)dZ (3.77) 

m J R 2jv 

where is the two-body density of /x. We write, for some parameter (3 to be fixed later on, 

^ + ^_2^ log|i _ 9|£(l _ w /M ! + y !VtogH _ logM _ l+rs {3JS) 



3 We do not change the notation for the scaled quantities. 
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with 

r,- U 'J* + .o 6 k-»l |. (3-79) 

Minimizing 

^(sc, y) := |a;| 2 + |?/| 2 - (5 log |x - y| 



with respect to y we find that the minimum is attained for \y\ = |- I 1 — y/l + 25/\x\ 2 ) and that 



/ ^ 3 1 1 2 5 1. . /p-fs — — \x\ + J\x\ 2 + 26 
<p{x,y) > ^N 2 - g - 2NVM 2 + 2 *-< yi °g — — *Y 



> \x\ 2 - £log I '"' ' v '"' ' ~ I > 



\x\ + y/\x\ 2 +2S 



\x\ 2 -6log(y\^T26) . 



Minimizing the last expression with respect to x we obtain 
from which we deduce, taking b = 4(N — l)/m/3, 



We can now choose 

N - 1 

= 2 

m 

and deduce 

r>>o. 

We can thus bound from below the free energy (|3.74[) as 

Fn[iA > FnM ■= [ H N {Z) f i(Z)dZ + — [ n(Z) log [i(Z)dZ (3.80) 

where 

N 

H N {Z) = ((1 - P) M 2 - 2 log |^| - 1) (3.81) 

3=1 

is now a one-body operator. The minimum i*V of J~n over symmetric probability measures fi is 
attained at 

/z = (g th )® N (3.82) 

with 

£ th (z) = -L cxp (-m ((1 - /3) |z| 2 - 2 log |z| - l)) (3.83) 

and 

F N = --\og(Z th ). (3.84) 

TO 



Moreover we have for any symmetric probability measure fx 

1 f 1 ^ 



F N [»]>F N + -I MlogT^r (3.85) 



>F N + - f M (1) log^, (3.86) 
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using the subadditivity of the entropy (see e.g. [Kiel Proposition 1]) for the second inequality. The 
equivalent of formula (|3.23p in rescaled coordinates reads 



1 



1 



e ( r ) = ^h ex P(--( r -21ogr-l) ), £ th = --logZ 



1 



7th 



(3.87) 



with Z th a normalization constant (we again keep the same notation for the quantities after the 
scaling of distances). In the new coordinates, g th resembles a Gaussian centered on the minimum 
of W m at r = 1 with characteristic length mT 1 ! 2 . Its maximum is thus of order m 1 / 2 and it decays 
exponentially fast in the region where |r — 1| 3> mT 1 ! 2 . Since j3 = 0(N/m) it is not difficult to 
realize, using p.83[) and (|3.87[) . that 

logZ th = logZ th + 

and thus it follows from (|3~M)l and (j3~5rJ)) that 

F N >NE th + — f /i (1) log^- + 0(N 2 m-' 2 ). (3.88) 
m Jr2 g th 

We now use (g th )® N as a trial state for and obtain 



* , 2iV(iV-l )j(gth)gth)- 



F N < NE 



Arguing as in the proof of Theorem 13. 11 Step 4, we approximate g th by a delta function along the 
circle of radius 1 to obtain 



g th h„ th = /i„ th (l) + 0(m- 1 / 2 ) 



But, using Newton's theorem (|3.30[) and the exponential decay of g th for \r — 1| 3> m 1/>2 , 
h gth (l) = -2n [ g th (r) (log r)rdr = O^ 1 / 2 ), 

Jr>l 

since the integral is located in a region where logr = 0(m -1 / 2 ) and g th is normalized. We thus 
have the upper bound 

F N < NE th + CN 2 m.- 3/2 , (3.89) 
which, combined with (|3.88[) and the CKP inequality p.27|) gives 

< CN^m- 1 / 4 



TV 



and there only remains to note that also 



„th ~th 



r"|lTV ^ CN l / 2 m^ A 



and scale variables back to deduce the desired result. 



□ 



Proof of (|3.18[) . In the course of the proof above we have established (still in rescaled variables 

N 

tf;v(Z)>^[(l-/3)N 2 -21ogM] 

3=1 
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and 

F N = — logZ N = NE th + 0{N 2 m~ 3 ' 2 ) 
m 

where Zn is the normalization constant of hm- In view of the expression (|3.83[) of g we have 



- log g th (z) = E th - (1 - 0) \z\ 2 + 2 log \z\ + 1 
m 



and thus 



1 N 
H N (zi, . . • , z N ) = — cxp (-mH N (Z)) < TT g th (z 3 )exp(N 2 m^ 2 ). 
Zn r~ 

3=1 

After integration over N — 1 variables and inspection of the expression for g th it follows that 

(Jff (z) < exp(— cm(r — l) 2 ) 
when |r - 1| > CNm~ 3 / 4 for C large enough. Then (|3.18[) is obtained by a change of scales. □ 

4 Improved energy bounds for the LLL problem 

In this section we use the results of Scction[3]to compute the energy of our trial states. This allows 
to improve the upper bounds in Theorem 12.31 More precisely, we prove 



Proposition 4.1 (Energy upper bounds). 

In the limit u, k — > 0, N — > oo we have 



ifuj> -2kN, 



E L < uN 2 {l + o(l)) + -fciV 3 (l + o(l)) (4.1) 
3 



£ L < -— + -kN 3 (1 + o(l)) (4.2) 



ifw < -2kN and \cj\ < kN 7 / 5 log N, and finally 



£ L <-^-j^V(l + (l)) (4.3) 

ifuj < -2kN and |cj| > ZsiV 10 / 3 . 

Proof. We use the trial states (|1.11[) . Since they all have zero interaction energy, we are left with 
estimating (remember the scaling of space variables in (]3.1[) ) 

r N 

Jm? N ]-[ 

= N ' 2 {^j (ur 2 + kNr 4 ) fj$(zUdz. (4.4) 

Our main tools are (|3.15[) and (|3.17p which essentially say that (j.^' w g cl or m g th , depending 
on the regime. The main terms in the right-hand sides of (|4.1[) and (|4.2p are obtained by replacing 
directly [iff by g cl or g th and our main task is to estimate the error. We cannot use (|3.15[) or (|3.17l) 
directly however because the norms appearing in the right-hand sides are certainly not finite for 
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V = bjr 2 + kNr 4 . We therefore first employ (|3.16p (or (|3 . 1 8[) ) to restrict the integration domain : 
We will use two smooth radial cut-offs functions Xin an d Xout satisfying 

Xin + Xout = 1 (4.5) 

and decompose £ L [ v f , „] as 

£ L [vK*] = N 2 f Xin V^ k (z)g(z)dz + N 2 f X inV u N k (z) (v$\z) - g(z)) dz 

+ N 2 f X outV u N k (z)^\z)dz (4.6) 

JR 2 

where g = g cl or g th depending on the regime and we denote 

V£ k (r)=u;r 2 + kNr 4 . (4.7) 

In the electrostatic regime our choice of cut-offs functions will ensure Xin = 1 on supp(g el ) and the 
first term is thus readily computed using the explicit expressions (|3.20[) and (|3.21j) : 

N 2 £ Xin {lot 2 + kNr 4 ) g°\z)dz = uN 2 (l + ^) + kN 3 (t + 2^ + ^ . (4.8) 



Optimizing the above expression with respect to m we find 



'-opt 



if to > -2kN 

-N - §: if uj < -2kN. 



(4.9) 



Therefore the Laughlin state is favored for to > —2kN whereas it is better to add a vortex at the 
origin for lu < —2kN. The term (|4.8[) becomes 



N 2 / Xin (ujr 2 + kNr 4 ) g e \z)dz = -N— + -kN 3 . (4.10) 
' — 4fc 3 



if ui > -2kN, and 



2 i 

N 2 I Xin (ur 2 + kNr 4 ) g e \z)dz = -N— + -kN 3 (4.11) 
J R 2 4k 3 

if lu > —2kN. In the thermal regime we keep the same expression for our choice of m and obtain, 
using (EE25D , 

f uj 2 3 

N 2 Xin (lot 2 + kNr 4 ) g th (z)dz < —N— — —ujN(1 + o(l)). (4.12) 
J R2 4k 2 

Note that the main terms above may be recovered from (|3.23[) by neglecting terms beyond quadratic 
in a Taylor expansion of W m around r opt . 

When estimating the remainder terms in (|4.6[) we distinguish between three regimes. The 
"cases" below refer to the different cases in Theorem [ 



Case 1 and 2, m opt < CN . We take g = g in (|4.6[) . The support of g is uniformly bounded in 
this regime. Using ()3.16[) we see that for r > i?+ 

M^W < Ccxp (-CN(r 2 - Clog AO) (4.13) 

and thus /z^ is exponentially small, both as a function of r and N for r > C(log N) 1 / 2 . We choose 

Xin = 1 in B(0,C(logA0 1/2 ), X in - out of B(0, 2C{\ogN) 1 / 2 ) 
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and x ou t accordingly. Then Xin = 1 on the support of g el as desired and (|4.8[) yields the main 
terms in the right-hand sides of (|4.1[) and (|4.2p . We also assume |Vxin| < C^logTV)^ 1 / 2 . We then 
note that 

|V(xinV5 fc )| < CilogN) 1 / 2 ^ + kN (log N) 3 ' 2 



for r < C(logTV) 1 / 2 while 



IXinWj^H^ < CQogN)\u\ + fcTV(logTV) 5 



Thus, using (|3. 15[) to estimate the second term in (j4.6[) we obtain that it is bounded above by 

C (loN 3/2 (log TV) 3/2 + kN 5/2 (log TV) 5/2 



and im follows, using (|4.10[) or (|4.11j) for the main term of P~o) . The third term in (jH))) is 
negligible thanks to (]4. 13[) . 

Case 3, TV < m opt < TV 7 / 5 log TV. We take again g = g cl in (|4.6[) . In this case (|3.16p gives 

M^M < Ccxp (-CN{{r - r opt ) 2 - Clog TV)) 
for r > i?+ or r < i?~ and thus [iff is exponentially small for \r — r opt | > C(log TV) 1 / 2 . Wc take 
Xm = 1 for r opt - C(logTV) 1 / 2 < r < r opt + C(logTV) 1 / 2 , Xin = for |r - r opt | > 2C(logTV) 1 / 2 

and Xout = 1 — Xin accordingly. As before the main terms of (|4.2[) come from (|4.8|) and the third 
term in (|4.6[) can be neglected due to the exponential decay of fi^K For |r — r opt | < C(logTV) 1 / 2 
one can easily realize that 

IVtenOI^CQogAO^M, 



whereas 



XinV^|| L2 < c/v^/^iog/v) 3 / 4 ^! 5 /^- 1 / 
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Using (|3.15|) . one can see that the error due to the second term in (|4.6[) is bounded by 

C/V 5 / 4 (log/V) 5 / 4 |w| 5/4 /€- 1/4 + CM/V^log/V) 1 / 2 
and this is negligible in front of kN 3 if |cj| <C kN 7 ' 5 log TV, which concludes the proof. 

Case 4, m opt ^> TV 10 / 3 . Here we take g = g th in (|4.6[) and notice that (|3.18l) ensures that fj,^ is 
exponentially small when 

k-r opt | > L := max(TV- 1/2 ,TV 1/2 m- 1/4 ). 

We therefore choose 

1 if \r - r opt | < L 
Oif \r-r opt \ >2L. 



Xin(r) = 



With this choice and using (|3.18|) . we can neglect the third term in (|4.6p since it will be much 
smaller than the two others. The main term is computed using (|4. 12[) . For the second term of 
gllD we use (|3~T7|) to see that 



N 2 I XiaV£ k {z)L%\z)-Q{z) 
;u>2 \ 



dz 



4.3 



where we approximate V^ k by its second variation around its minimum. The last quantity is 

0(JV 7 / 2 A; 3 / 4 M 1/4 ) 
when \io\/k oc m < TV 4 (recall the choice of to = TO opt in (|4.9[0 and 

0(7V 3 / 2 fc 1 /4 w 3/4 ) 

otherwise. We conclude that the error term is negligible in front of the subleading term —^ojN in 
(14. 3[) when \u>\ ^ fcY 10 / 3 , which concludes the proof. 

□ 

For completeness we reproduce the argument from [RSY] allowing to see that at least the 
order of magnitude of our energy upper bounds is correct. This is rather simple: We consider the 
problem of minimizing the potential energy in a sector of given angular momentum. This makes 
sense because we are in the lowest Landau level and we can rewrite the potential energy as in 
(|2.13[) . Let us denote 

E (L) = inf J (f, J2 h 3 F ^j > c nF = LF, (F, F) BN = 1 1 (4.14) 



where 



h = (w + 3k)zd z + k{zd z ) 2 (4.15) 



and hj is the same operator acting on the variable Zj. 
We have the following lemma : 

Lemma 4.1 (Potential energy at given momentum). 

For any L £ N 

L 2 

E (L)>(Lo + 3k)L + k—:=e(L) (4.16) 
with equality if L is a multiple of N . 

Proof. Since Cn and commute, they can be diagonalized simultaneously and we have the 

operator inequality 

' N 




from which (|4.16[) follows. The fact that there is equality when L is a multiple of N is proved by 
taking the trial state (tt^!)-^ 2 nj=i Z f /N - □ 



To see that the lower bounds of Theorem 12.31 follow, simply use this lemma with L = Lq, the 
momentum of a ground state, as estimated in Theorem 12. II 
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